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Abstract 



We derive the asymptotic behavior of determinants of truncated Wiener-Hopf 
operators generated by symbols having Fisher-Hartwig singularities. This task is 
achieved thanks to an asymptotic resolution of the Riemann-Hilbert problem associ- 
ated to some generalized sine kernel. As a byproduct, we give yet another derivation 
of the asymptotic behavior of Toeplitz determinants having Fisher-Hartwig singulari- 
ties. The Riemann-Hilbert problem approach to these asymptotics yields a systematic 
although quickly cumbersome way to compute their sub- leading asymptotics. 
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1 Introduction 



There is a long history of understanding the asymptotic behavior of determinants 
of structm-ed matrices. It started in 1915 with the seminal work of Szego on the 
asymptotic behavior of large size Toeplitz matrices. His result is known today as 
the strong Szego theorem. It states that for non-vanishing and regular functions 
b€C^ (T,M+) 

/c]0 



(1.1) 

The constants G [b] and E [b] are expressed in terms of the Fourier coefficients of log 6: 

G[6] = e['-^]= E[b]=eS.'''''''^'"'''-^ [log6], = [ ^e^^^ lo,b (6) . (1.2) 

J 27r 



The result of Szego underwent several refinements. In particular Baxter [6], Hirschman 
[26] and, finally, Ibragimov [27] successively weakened Szego's original assumptions 
on b. Furthermore, Widom |i39j considered determinants of block Toeplitz matrices 
and provided a clear interpretation of the constant E [b] in terms of an operator 
determinant. Despite the possibility to consider matrix valued functions b in the Szego 
theorem, there are limitations of its applicability, even in the scalar case. Indeed, 
the theorem already breaks down in the case of symbols having zeros, power-law 
singularities or even jump discontinuities on the unit circle Such symbols can be 
represented as 



a {6) = b{e)Y[us^^^^{e''/ap) with G ^ , (1.3) 
p=i 

^i{9—nsgnd)Sp 

(2~- 2 cos ( 



and u:s,rf,U') = o.^.m> ,0e]-7T-7T[ . (1.4) 



In such a representation one assumes that the function b is regular enough, non vanish- 
ing on and has a vanishing winding number. The conjecture about the asymptotic 
behavior of Toeplitz matrices generated by such symbols goes back to Fisher and 
Hartwig in 1968 [23| . More precisely, they claimed that 

n 

det^[T{a)] ~ (G[&])'"m£ C ([6], m^{7K) (1.5) 

m— s>+oo 

The Fisher-Hartwig conjecture was indorsed by a couple of examples where the authors 
were able to compute the Toeplitz determinants explicitly. The value of the constant 
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was later conjectured to be equal to 



n 



cm,{5ri,m) = E[b]ii 



Gil-j,-6i)Gil-j, + 5i 



i=l 



G(l-27, 



xUb'H-'^ ia,)bl^+'^ {a,)l[{l - ap/a.f"^^-^^'"-^^'^ . (1.6) 

i=l p^q 

There G is the Barnes' function and b± are the Wiener-Hopf factors of b, ie b = 
b+G [b] b- with 6+, resp. 6_, a holomorphic function on the interior, resp. exterior, of 
the unit disk and such that 6+ {z = 0) = 1, resp. 6_ (z) — > 1. 

The above conjecture was first proved for some particular cases of the parameters 
and 17^. Basor [2j and, independently, Bottcher [7] treated the case of several 
jump discontinuities (Vp ,7p = 0) under the restriction |3?(5p)| < 1/2. In 1985, 
Bottcher and Silbermann flO] proved the conjecture in the case |K((5p)| < 1/2 and 
|3?(7p)| < 1/2. Finally, Ehrhardt and Silbermann [21j proved the conjecture in the 
case of a single Fisher-Hartwig type singularity for all ranges of parameters 6p and 
7p where it made sense. This allowed to prove the conjecture in most of the cases 
involving multiple Fisher-Hartwig singularities [20J . The proof was based on the so- 
called separation technique developed by Basor [Sj. The Fisher-Hartwig conjecture 
breaks down in the case of the so-called ambiguous symbols. Basor and Tracy [3] raised 
a generalized Fisher-Hartwig conjecture for the behavior of Toeplitz determinants 
generated by such ambiguous symbols. This conjecture was proven recently in the 
framework of Riemann-Hilbert problems for orthogonal polynomials by Deift, Its and 
Krasovsky [TilfTS]. 



There exists a continuous analog of Toeplitz determinants, the Fredholm determi- 
nant of truncated Wiener-Hopf operators. The underlying operators act on functions 
g G (M) according to the formula 

X 

(1 + K) .g (t) =g{t) + j dt' K{t- t') g {t') dt' . (1.7) 



The kernel K is traditionally defined in terms of its Fourier transform a—l,ieK{t) = 
[<T - 1] (t). In the foUowing, we choose the below convention for the restriction of 
the Fourier transform (and of its inverse) to (M) nL^ (M): 

b] (t) = iO e~^*« and T [h] {() = J dth (t) e**« . (1.8) 

So that, for functions (a — 1) € (M), we have an explicit integral representation 

Kit) = ^ /dek(e)-l]e-^*f . (1.9) 
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The question of the x — >■ +00 asymptotics of det [/ + K] were first addressed by 
Achiezer [1] and Kac |31) . More precisely, they showed that for symbols a regular 
enough 

det [I + K] ~ exp < 



hoo 



^00 

log [a di+ j i log [a {i)] log [a (-^l d^ 



(1.10) 

There exist many generalizations of this formula. These either extend the result to 
less regular or matrix valued symbols a. However, just as in the Toeplitz case, the 
theorem breaks down when a has some jump discontinuities or power-law behavior. 
The continuous analogue cr of a symbol with Fisher-Hartwig type singularities reads 

-{i) = F{i)^a^,,.,{i-a,) , -.-iO=[^^J[^-^J . (1.11) 

Note that the exponents in the definition of cTu^u should be understood in the sense of 
the principal branch of the logarithm, ie arg € ] — vr ; vr [. In the above decomposition 
for (7, the function F is supposed regular enough and CTuf,,v^, (0 h^-s a singularity at 0: 

gi7r<5fcSgn(5) 

'^i^fe.i^fe (6 ~ — TTTo^T — with 27fc = i/fc+I7fc and 26k =Vk - i^k ■ (1-12) 

There exists a continuous analogue of the Fisher-Hartwig conjecture for Toeplitz de- 
terminants and it is due to Bottcher [8]. It was inspired by the study of truncated 
Wiener-Hopf operators that are generated by rational symbols; indeed, in the latter 
case, the author was able to estimate the Fredholm determinants explicitly. This 
conjecture was confirmed in many particular cases: the case |12j of pure jump type 
singularities, the case [9] where all = or all = and finally in the case [4] 
of a single pure Fisher-Hartwig singularity a^^v under the restriction (z^ it I7)| < 1. 
All these results were established thanks to some identity relating determinants of 
truncated Wiener-Hopf operators to determinants of Toeplitz matrices, and then the 
use of the formulae for the asymptotic behavior of Toeplitz determinants. 

Yet there exists an alternative approach to the asymptotic analysis of large size 
determinants of structured matrices. In true, it is well known that a Hankel matrix 
can be expressed as a product of the leading coefficients of the orthogonal polynomials 
with respect to the measure defining its entries. It was observed by Fokas, Its and 
Kitaev [24j that one can recast the problem of computing orthogonal polynomials into 
a certain Riemann-Hilbert Problem (RHP). This problem can be solved asymptoti- 
cally for a large case of weights |17 |[T8l[36] . As noticed by Krasovsky [35) . one can 
relate the asymptotic solution of the RHP for orthogonal polynomials with respect to 
a weight having a finite number of power-like singularities to the asymptotic behavior 
of the Hankel determinant defined in terms of this weight. Its and Krasovsky used 
an analogous identity to establish the asymptotic behavior of a Hankel determinant 
defined in terms of a gaussian weight having a jump discontinuity [29]. Moreover, 
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still in the framework of RHP for orthogonal polynomials, Krasovsky estimated the 
asymptotic behavior of Toeplitz matrices on an arc and generated by symbols having 
jump type discontinuities [34J thanks to the relationship between polynomials orthog- 
onal on an arc and those on a line segment. His approach presented no obstruction 
for a generalization to the case of root type singularities. The case of Toeplitz, Hankel 
and Hankel+Toeplitz determinants with Fisher-Hartwig singularities has been treated 
recently in the framework of the Riemann-Hilbert problem for orthogonal polynomials 
by Deift, Its and Krasovsky |14yi5) 

Truncated Wiener-Hopf operators being continuous analogs of Toeplitz matrices, 
there arises a natural question concerning a RHP approach to study the large x be- 
havior of Fredholm determinants for such operators. In this paper we will show how to 
tackle, in the framework of RHP, the x +00 asymptotics of determinants of trun- 
cated Wiener-Hopf operators generated by Fisher-Hartwig symbols. This treatment 
is based on a relationship between truncated Wiener-Hopf operators and the so-called 
generalized sine kernel acting on M. The latter kernel is an integrable integral operator. 
Such operator can be analyzed by a RHP as observed in |3^. We asymptotically solve 
this RHP. The construction of its asymptotic solution is an extension of the work |33) , 
the latter being itself a generalization of an unpublished study on the pure sine kernel 
by Deift, Its and Zhou. This approximate resolvent allows to compute the leading 
asymptotics of the Fredholm determinants of truncated Wiener-Hopf operators. The 
latter constitutes the main result of this article: 

Theorem 1.1 Let I+K he a truncated Wiener-Hopf operator acting on the segment 
[0;x] and generated by the symbol (a — 1) £ (M) with 

n 

c^iO=F{0'[l(r^„v,{C-ap) ateR, ai < • • • < a„ (1.13) 

k=l 

where 

• F is non-vanishing and holomorphic in some open neighborhood U of the real 
axis ; 

• F — 1 G L (M) and moreover F {^) — 1 = O 2 j ^ for some k > and 
^ +00 in U ; 

• |K(5fc)| < 1/2, andK(7fc) < 1/4 /or A: € [ 1 ; nj. 

Then the leading asymptotics of the 2-regularized determinant of I K read: 

n 

-pr Fl" (ofc) F^" (ofc) -pr / (ofc - gp if 

M K' (flfc + ^- ("'^ - ^) l}p V ("'^ - «p + ("'^ - "p) 

(1.14) 



'1 + 



")) 
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There 



Gl [a] = exp L I ^ [log [a (0] + 1 - ^ (61 [ , (LIS) 



E [F] = exp < 



+00 I 

j d^CJ"-' [logF] (0 [logF] (-0 \ , (1.16) 

J 



and F± are the Wiener-Hopf factors of F : 

» 

T/ie estimates for the correction involve the constant p = 2max|5R ((5fc)| < 1. Finally G 
stands for the Barnes G function and we remind that det2 [I + K] = det [(/ + K) e~^] . 

Note that a similar result can also be established for symbols a — 1 £ (M). 
The above theorem reproduces all of the aforementioned results that were obtained 
for particular cases of singularities. However, it shows that the original conjecture 
doesn't hold in its whole generality. Indeed the conjecture [SlfTT] predicts the presence 
of 

[ok - apY + 1 



/ (n, - n ^^ 



n 

k<p 



:i.l8) 



whereas we find 




:i.i9) 



Of course both results coincide in all the cases previously investigated ie = ii^fc, 
V/c , ffc = or VA; = 0. Moreover our approach opens a way, at least in principle, 
to asymptotically inverting truncated Wiener-Hopf operators generated by Fisher- 
Hartwig symbols. Such an inversion could be carried out in the spirit of the inversion 
for holomorphic symbols a proposed in |33| . 

As a byproduct, using a formula [16] relating certain Fredholm determinants of 
a generalized sine kernel to Toeplitz matrices, we derive the asymptotic behavior of 
Toeplitz determinants generated by symbols having Fisher-Hartwig singularities. This 
computation reproduces the result of Ehrhardt [20j for the leading asymptotics and 
of Deift, Its and Krasovsky for the sub-leading ones |l5j. Hence, we see that Toeplitz 
determinants and those of truncated Wiener-Hopf operators are both related to a 
determinant of a generalized sine kernel. The only difference being the interval on 
which the generalized sine kernel acts. 

The article is organized as follows. In the first Section we establish the link between 
Fredholm determinants of truncated Wiener-Hopf operators and those of a generalized 
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sine kernel. The second part of Section [2] is devoted to introducing some notations. 
In particular, we shall consider two types of symbols for K. The first one belongs to 

(M) whereas the second to (M). Both have the same Fisher-Hartwig singularities, 
but in the second case, part of the function F (jl.lip contains some prefactor depending 
on 5 in order to ensure that (u — 1) € (M). 

In Section [3] and S] we asymptotically solve the Riemann-Hilbert problem associ- 
ated with the generalized sine kernel. In Section [5] we derive the asymptotic behavior 
of the resolvent of the generalized sine kernel. In Section [6] we use this asymptotic re- 
solvent to obtain the asymptotic behavior of the Fredholm determinants of truncated 
Wiener-Hopf operators under investigation. We compare our results with the already 
existing ones. In the last Section we adapt the RHP associated to the generalized sine 
kernel in order to study the asymptotic behavior of Toeplitz matrices with Fisher- 
Hartwig symbols. We obtain the leading asymptotics of such Toeplitz determinants 
and also compute their first sub-leading corrections. The structure of these corrections 
indicates that at least part of the asymptotic series can be deduced by making shifts 
of certain parameters appearing in the leading asymptotics. This leads us to raise a 
generalization of the Basor-Tracy conjecture. 

2 Some preliminary results and definitions. 

In this Section we establish the link between truncated Wiener-Hopf operators 
and generalized sine kernels in the case of general and symbols a — 1. We also 
present some determinant identities that will be useful in our proofs. 

2.1 Truncated Wiener-Hopf and the modified sine kernel 

Let I -\- K he a truncated Wiener-Hopf operator acting on (M), ie. I + K acts 
on functions g E (M) as follows 

X 

(I + K) .g (t) = g{t) + j dt'K {t - t') g {t') dt' (2.1) 



It is useful to define the kernel K in terms of its Fourier transform cr — 1, ie K (t) = 
T-^ \a-\\ (t). 

We shall focus on two cases of interest: a {^) — 1 E (M) and cr (^) — 1 G (M). 
In the first case, I + K is trace class and hence its determinant is well defined. In the 
second one, I -\- K is Hilbert-Schmidt; hence, one ought to consider the 2-regularized 
determinant ^oil + K ie det2 [I + K] = det [(/ + K) e"^] . For the purpose of this 
section only, we introduce the index I in (T(^). This means that (c^^) — l) € (M), ie 
(cT(i) - 1) G (M) and (a(2) - l) G (M). 

Lemma 1 Let K(^£-^ = [^{e) ~ l] > then one has the identity 

I + K^i) = F'^ oMoDo[l + V(i)) o D-^ o Ar^ o T . (2.2) 
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Similarly, 

The operators I + V(i), resp. I + V(2), act on (M) with kernels 



27ri — r/) 



6 2 _ e 2 



6 — 17 



2i7r — rj) 

and D, M are the multiplication operator on (M) 

{D.g){i) = r-^ g{i) . 



Proof — In any of these two cases, we have that 



(2.3) 

(2.4) 
(2.5) 

(2.6) 



J" 



K^p) {t - t') g {t') dt' 



(0 



2ilT (^ — 7]) 



[OpoV^p^oOp' oT][g]iO . (2.7) 



Where Oi = M o D and O2 = D. Here, we precise that Vj-i) o O^^ is indeed a weh 
defined operator on (R) . □ 



These two identities relate the truncated Wiener-Hopf operator to a generahzed 
sine kernel acting on M. If one is able to construct the resolvent for this operator, then 
one is able to invert the corresponding Wiener-Hopf operator by using (12. 2p or (12. 3p . 
This correspondence has already been used in [33] to build the resolvent of truncated 
Wiener-Hopf operators whose symbols are holomorphic, non-vanishing functions on 
some strip around the real axis that are decaying fast enough at infinity. In the case 
of symbols ai having Fisher-Hartwig singularities as in (II. lip , the expression for the 
leading asymptotic resolvent of the underlying generalized sine kernel is much more 
involved than for operators considered in ^33j, hence taking the Fourier transform and 
then obtaining some manageable result might be complicated. 

The two identities given in Lemma [2] allow to establish a connection between suffi- 
ciently regularized Fredholm determinants of / + K(^£-^ and those of the corresponding 
generalized sine kernel. Hence, to study the asymptotics of det^ [/ + it is enough 

to focus on the ones of the associated generalized sine kernels. 



Lemma 2 Let Kt 
(I2II), then 



in 



ll be the kernel of the integral operator given in 



det2 [I + K(2)] = det2 [I + T/(2)] , 



(2.8) 



8 



and 

det [I + = det [I + . (2.9) 

Proof — The second equality can be obtained thanks to the Fredholm's series repre- 
sentation for the determinant of a trace class integral operator: 

hoo 



det [/ + = ^ / d"t det„ {U - tj)] 

n=o i 

+00 -II- p rl"',^ " 

= E / / ^ n (-(1) (^P) - 1) 



n=0 n K P=l 



+CX3 „ 

5^- / d"edet„[y(i)(Cp,e,)] (2.10) 

n=0 



Now, let us prove the first identity. Define i?2 (^(2)) ~ (-^ + -^(2 

)) e-^(2) _ J, Since 

K2 is Hilbert-Schmidt, we have that R2 (^{2)) is trace class |37j. Moreover, equation 
(|2.3p implies that 

(K(2)) = 0D0R2 (^(2)) o D o J- (2.11) 

and i?2 (^(2)) is trace class as V(2) is Hilbert-Schmidt. Moreover, the Fourier tansform 
J-" its inverse J-""^ as well as D and Z)"^ being continuous operators on (M) we have 
that R (V(2)) o -D^^ o is trace class. We can thus change the order in the operator 
product appearing in the determinant so that 

dei2[l + K(2)\ = dei[l + T-^oDoR2{Vi2))oD-^oT] 

= det [/ + i22 (V(2)) oD-i oj-o oD] = det [/ + i?2 (%)] □• 

2.2 General assumptions 

Motivated by the latter results, we consider the two Fredliolm operators I + Vi 
acting on (M) and defined through equations (12. 4p and (12. 5p . The operators are 
defined in terms of the symbols cr^ below 

n 

<y{() iO = Fn) iO n ^-^'^^ - "'^) • (2.12) 

k=l 

The functions Fr£\ (^) are chosen according to 



FmiO = &(i)(e)n(i + ^) , (2.13) 

Fi2)iO = &(2)(0 • (2.14) 



and we assume that the functions are such that 
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• is holomorphic on some open neighborhood U of M; 

• never vanishes on U ; 

( l + K \ 
j^l ^ ] for some k > 0, when ^ — ?> oo in [/. 

We also make some assumptions on the exponents 25i =Vi — vi and =Vi + uf. 

• \/k, K(7fc) < 1/2 in the (M) case; 

• V/c, 3?(7fc) < 1/4 in the (j^) case; 

• V A; , |K < 1/2 independently of the (M) or the (R) case. 

The behavior of a^^u (0 around = shows that the last restriction on 6k covers 
almost all the possible types of jump singularities (T(£) could have. However, the 
case = ±1/2 for some k's should be treated separately. In particular, one 

expects additional corrections to the asymptotic formula ()6.36p . These should have 
the same structure as those appearing in the Basor- Tracy conjecture for Toeplitz 
matrices [3]. We have the 

Lemma 3 Under the above assumptions, the symbols cr(f) given in ()2.12p are such 
that o-(i) - 1 G (M) and a^2) - 1 G (M). 

Proof — We give the proof in the (M) case only. The assumptions on the param- 
eters 7fc and the local behavior of cTu^;u^ H — o-k) together with oi < ■ ■ ■ < enure 
that (T(i) — 1 G {[~^ ) -^]) foi' any finite M. It remains to check the integrability 
at infinity. (yuy,Mk decreases at infinity as: 

^u,,v, (0 = 1- 2z4r ' + O (r') (2.15) 

Therefore, we have 

( . n \ " 

^ k=l J k=l 

= 6(i)(0-i + o(r') = o(iei'+'^) , (2.16) 

and the claim follows. □ 
2.3 The resolvent 

Let /jf be the solutions to the integral equations 

/« (0 + / (,) - 1) f^;||^/(.. (,„ d,, = e-i . (2,17) 

R 
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It is well known [30] that the resolvent operator / — ^(i), resp. / — ^(2), of / + V(i), 
resp. / + V(2), has a simple expression in terms of Indeed 



Y^c^(i) (0 - Iy^o-(i) iv) - 1 



2i7r — ij) 



/f iOf^y {v)-r+' iv)/'-' (0(2,18) 



(1) („\ 



^(2) io - 1 

2i7r — T]) 



(2) /^A f (2) 



(2.19) 



Lemma 4 Suppose that det [/ + V(i)] 7^ 0, resp. det2 [l + V(2)] 7^ 0, i/ien i/ie solu- 
tions of (I2.17P are entire functions. 

Proof — Suppose that det [/ + V(i)] 7^ 0. Then / + V(i) : L'^ (M) (R) is 

invertible and its inverse / — : (R) — ?> (M) can be constructed in terms 

of a Fredholm series. Thus, since e^*^2 \j ('^) ~ ■'■ ^ the unique solution 

7± = ./cT(i) (6 - 1 (6 of 



/± (e) + j Vn) iC, V) f± iv) dT] = e±*4 ^a(i) (0 - 1 (2.20) 

belongs to L2 (M). We also have that 

VeeC, ry^^cxd) • (2-21) 

Therefore, 

. r? ^ F(e,r?) ^ /£) (^) (^(^ (^) _ 1) G (M) , for ah ^ G C; 
• C i-T' (^, ??) is entire for almost all r] . 

Thus, j F{i,r])Ar] (2.22) 

R 

is an entire function. By (|2.17p . so is In the case det2 [/ + V2\ 7^ 0, we have 

immediately that (l — (T{2) (0) /± ^ ^ (^)' the rest of the proof goes the same. 
□ 

2.4 Determinant identity 

We shall now derive an important determinant identity. This identity allows to 
obtain the leading asymptotics of truncated Wiener-Hopf operators from thoses of its 
resolvent. 
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Lemma 5 Let V^i^ and V^2) ()2.4p and (12. Sh and defined in terms of the 

symbol (I2.12p subject to the assumptions of subsection \2.^ We also assume that 
det [/ + V(i)] 7^ and det2 [/ + ^(2)] 7^ 0- Suppose that (5p equals 6p or jp, then the 
following identities hold 

d^^log det [I + V^,)] = [ ^^'][f^\ dp^a(^i){Od^ , (2.23) 

J o-(i) l?j - 

R 

d^Jogdet2[l + V^2)] = ||^^^||^a^,a(2)(6-5;3,%(e,0}^ (2.24) 

Proof — We first treat the L (M) case. Let {5p,7p}"^^ be a point in C^" fulfilling 

the assumptions of subsection (12. 2p for the (M) case and such that det [/ + V(i)] 7^ 
0. It then follows from the Fredholm series for det [/ + V(i)] , that the latter is a 
holomorphic and non-vanishing function on some open neighborhood of 7p}p_-^. 
It is in particular differentiable and its derivatives can be expressed by using the 
resolvent operator / — R(i)- Setting e'^'^^^^^ = (7(x) {(,) — 1, we get 

dp^ log det [/ + 1/(1)] = tr { (/ - i?(i)) . dfs^V^,^ } 

= Jdc dp^Gi (e) [^(1) • (/ - i?(i))] + [(/ - R(i)) ■ V(i)] (e,e) dp^c, (o 



2 j de {u) dp^G, io = /de i?(i) {to ^^^§^ 



In the intermediary equalities we used the symmetry of the kernels as well as the 
cyclicity of the trace. The (M) case is proved by density. Let Xe be the characteristic 
function of ] -e ; e [. Then I + V2;e, with V2-e v) = Xe (0 V(2) iC, v) Xe (v), is trace 
class for all e > so that 

det2[I + V2;e] = det[I + V2-e]e-''^'-'^- . (2.25) 

As det2 [/ + ^(2)] 7^ 0, det2 [/ + V2-e] 7^ for e large enough, and hence det [/ + V2-e\ 7^ 
as well. One can then apply the results for (R) kernels for the /3p derivative. We 
get, 

dp^ log det2 [/ + V2;e] = JdXxe (A) ^^^^^^dp^a2 (A) - J dAxe (A) dp^V^2) (A, A) . 

M R 

(2.26) 

The e — >■ +00 limit in the rhs becomes licit after merging the two integrals into one. 
□. 

The asymptotic solution of the RHP presented in the upcoming sections will allow 
us to construct approximate in x resolvents of Vi and V2 uniformly in respect to the 
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parameters 5p and jp. It will then remain to use these approximations to compute, 
in the large x limit, the integrals appearing in (j2.24p or (|2.23p . Once this is done, 
it is enough to integrate the result from Pp = to Pp. Such an integration is analo- 
gous to the separation technique in the operator approach to asymptotics of Toeplitz 
determinants p]. It is then not a problem to repeatedly apply the procedure so as 
to obtain the asymptotics of the determinant. At this stage, it becomes clear why, 
in our approach, the (R) case doesn't follow from the (M) one. As a matter of 
fact, if we want to keep jump singularities and still have an (R) kernel, we ought 
to add an additional factor depending on the (5p's as it was explicitly done for the 
function cf (|2.12p . This modifies the 5p dependence of the integrand in (j2.23p 

and hence the integration procedure. The result should also be, in principle modified, 
but eventually we see that the (R) case can be obtained from the (R) one by 
restricting correctly the parameters and replacing i^(2) by ^(2)- However, since these 
are only minor modifications, from now on we focus on the (R) case. The interested 
reader will find no problem in adapting the proofs to the (R) case. Accordingly, 
from now on, we drop the i subscript labeling (T(£), the kernels V(£) and the resolvents 
R(^iy We will denote these quantities by a, V and R and assume that a = ct(2) 
defined in (l2T2]> . 

3 The Riemann Hilbert Problem 

We start this Section by introducing a RHP for a matrix x- This type of RHP 
is adapted for constructing resolvents of integrable integral operators |30| such as the 
generalized sine kernel. We then perform a few transformations of this original RHP, 
in order to boil it down to one where the jump matrices will be /2 + o (1) uniformly 
away from the points and in respect to the x +oo limit. The first step will 
consist in finding a scalar valued function a such that x^^'^^ has a jump matrix with 
1 in its lower diagonal entry. Then we deform the original cut. The jump matrix on 
the new contour has the desired properties. 

3.1 The initial Riemann-Hilbert problem 

As first observed in [30j, the problem of finding the resolvent of any integrable 
integral operator is equivalent a Riemann-Hilbert problem (RHP). Indeed, let f± be 
the solutions of (I2.17P for the corresponding a. Then, there exists [30j a matrix x 
allowing to reconstruct the solutions f± of the integral equation (|2.17p , and hence the 
resolvent: 




e-4,e«l)x-^(e) • 
(3.1) 



This matrix x solves the RHP: 



• X is analytic on C \ R 
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VA: € 1 1 ; nj , there exists G GL2 (C) such that 
X = Mk<l2 + g{z)Bk + \z -ak\{giz) + 1)0 (^^ | ]} , z ^ ak 



z^oo 



1 
1 



. x+iz)G{z) = X-{z) ; zGR . 

The matrix B (z) appearing in the estimates around is a rank one matrix that 
takes the precise form 



Biz)=[^ ^ J . (3.2) 

The function g reads 

, . f ds a (s) — 1 , . 

R 

and has the local behavior 1251 at z — ?• a^: 



SW= 0(l) + 0((.-a.r^-) tor ,,^0 
i O (log {z - Ofc)) for 7fc = 

The matrices are apriori unknown and will be determined once the solution is 
known. What only matters for the solvability is the invertibility of M^. Lastly, we 
adopt the convention that the symbol O (M) for some matrix M is to be understood 
entry wise ie x = O (M) means that Xij = O (Mij). 
The jump matrix G appearing in the RHP reads 

G(z)-( (a(z)-l)e-^\ 

(i-a(z))e-*"^ a{z) J ■ 

Finally, x+ (i) (resp. x~ (^)) stands for the non-tangential limits of x (z) as z ap- 
proaches a point t of the contour from its + (resp. — ) side. 



+ 



Figure 1: Original contour for the RHP. 

Proposition 3.1 Whenever ^(7^) < 1/2 , VA; € [1; and det2[I + K] / 0, the 
solution to the RHP exists and is unique. 
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Although most of the proof is standard, see for example |36j (careful handling of 
the singularities) and |13| (general exposure) , we include it for the sake of completeness 
and the reader's convenience. The last part of the proof dealing with a cancelation of 
the singularities due rank(i? (z)) = 1 is new. 

Proof — 

The RHP is equivalent to the singular integral equation for x- 

Since we assume that det2 [/ + -ftT] ^ 0, as already mentioned, the resolvent operator 
I — R exists. Hence, one can express a solution of equation (13. 6 p in terms of R or, 
equivalently, in terms of f± ()2.17p (whose existence follows from the existence of the 
resolvent): 

= W ^"SH^ ( Y> ^:::! l'\ ti ) ^ (3.7) 

J 2nT z - s \ -/_(s)e "2 /„(s)e"2 J 

R 

This proves the existence of solutions provided we show that (|3.7p has the desired 
behavior around each point a^. The solution of the jump conditions given in (13. 7p can 
be written as 

XW = + (-^!t-2 ^!t2)/^^^ (3.8) 

V -f-{z)e 2 f-{z)e"-2 J J 2nT z-s 

R 

with 

-<^)--/^^{(;!l:i)(-"*--" 



-e ^^2^e^^2 



The matrix M (z) is holomorphic around z = Uk, and it is easily seen from (|3.ip 
that 

Therefore, 

ds a (s) — 1 



x{z) = M{z)lh + g{z)(_i^^ \ with g{z) 



2in z — s 



(3.10) 

The local behavior of g around z = can be inferred from [25] . Then, the claim for 
the local estimates follows by expanding the holomorphic matrices around z = and 
setting M (a^) = M^. We shall now prove that Mk € GL2 (C) 
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Given any solution x to the above RHP, det [x] is analytic on C\M. Since det [G] = 
1 we have that det [x] is continuous across M \ U^^-^ {^p} ^^'^ can thus be extended to 
an analytic function on C\UpL^ {<^p}- Using the estimates for x and the multilinearity 
of the determinant, we get that for z ^ 

det [x] = det [Mfc] det [h + g {z) B (z)] + O {g {z) {I + g {z)) \z - a^l) . (3.11) 

As B {z) is a rank one matrix det [I2 + g{z) B {z)] = g {z) tr [B {z)] = 0, and thus 
the first term is a 0(1). The estimates for the local behavior of g{z) together 
with the hypothesis on the parameters 7^ ensure that O {g {z) {1 + g {z)) \z — a}~\) = 

o (^z — a^y^^. det [x] has thus no pole at z = a^. Its singularities at the a^'s are 
thus of a removable type and hence det [x] is holomorphic around a^. It follows that 
det [x] an entire function that is bounded at infinity in virtue of the normalization 
X — > l2- By Liouville's theorem, det [x] = 1. In particular det [x] (ofe) = 1, which 
can only happen if det [Mfc] = 1. 

We end the proof by showing the uniqueness of solutions. Let xi and X2 be two 
solution of the RHP for %. Then as det [X2] = 1; X2 is analytically invertible on 

n 

C\M. The matrix XiX2^ is holomorphic on C\M and continuous across ]R\ IJ {ak}- 

k=l 

Moreover, using the local behavior at z = we get 



X2' 



(I2 +* Comat (B (z)) g (z)) M^^^ +o{{l + g {z)) \z-ak\{^ \ \ 



(3.12) 

Here, we used the fact that the inverse of X2 is given by the transpose of its comatrix 

due to det [x2\ = 1- We have also introduced two, a priori distinct, matrices M^^^ 
associated with each of the solutions xt- Computing the matrix products and using 
that B {z) +* Comat {B {z)) = 0, B {z) •* Comat {B (z)) = 0, wc get 

Xi{z)X2H^) = Mk,iM^^l + o(g{z){l + g{z))\z-ak\(^ \ J , (3.13) 

The local estimates for g imply that g (z) {1 + g (z)) \z — ak\ =0 (^\z — a^p^ j . Hence 

Xi (^) X2 ^ (^) no poles at z = ak- The singularities at the a^'s are thus removable 
and, because of the asymptotic condition, we have XiX2^^ = I2 ■ This guarantees the 
uniqueness of the solution to the RHP, at least for ^{"fk) < 1/2- O 

Note that one doesn't have to make such fine estimates for proving the uniqueness 
of solutions if one would assume that ^{'jk) < 1/4:. However, we presented here this 
more complex proof as it also holds in the (M) case. 

3.2 A helpful scalar Riemann-Hilbert problem. 

Let a be the solution of the following scalar Riemann-Hilbert problem 

n 

a is analytic on C \ M , a_ (^) = a+ (^) a (0 , 4 € M \ I J {ofe} , a (0 — ^ 1 • 

(3.14) 
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This Riemann-Hilbert problem can be solved almost explicitly for the particular form 
(jl.lip of (T we study, namely, a (z) = (z) for z G and a (z) = (z) for z G Ti-, 
where 

at (^) = F^' (z) n ^ e , 

fc=l - fflfc + 

=F-{z) ft ( '""'"T " -e^- . 

fc=i V z-ak J 

'H± is the upper/lower half-plane and F± are the Wiener-Hopf factors of F, ie F = 
Fj^F- with Fj^ (resp. F-) analytic in the upper (resp. lower) half-plane and going to 
1 at z — 7> oo in T-Lj^ (resp. 71^). There is no constant factor Fq in the Wiener-Hopf 
decomposition of F as F — > 1. 

The Wiener-Hopf factors of F have an integral representation either in terms of 
Cauchy or Fourier transforms of log F: 

logF+(z) = J ^^lS^ = F[Z{OF-^[\ogF]{0]{z) ,zGn+, 

logF_(z) =-J^^L^^ = ^[Ei-OJ'-'[logF]{0]{z) ,zGn., 

and H is Heaviside's step function. As F is analytic and non zero in U, it is clear 
from these integral representations that F^ and F^ have an analytic continuation to 
U. Moreover, as F is non-vanishing in U and the decomposition F = Fj^F^ is still 
valid on [/, and F- have no zeroes on U . 
We introduce the auxiliary functions 



ap(z) = F{z)J[ 



k=l 




^ n U^^) (^^) e-M,(3.15) 

k=p+l 

dp{z) = {z - apf'^P ap{z) . (3.16) 

The functior0 Up, resp. dp, is holomorphic on {z : < 3fJz < ap+i}n[/, resp. -Dap,e = 
{z S C : |z — ttpl < e}. Here and in the following, e is such that -Dap,e C C/ and 
-Dap,e n F)aq,e = for p ^ q. The functions fip and CTp can be though of as the analytic 
parts of the local behavior of a on the real axis. Namely, one can continue a by 



^We stress that the subscript p appearing in ap has nothing to do with the notations of section 2. 
In the following Op will always refer to the definition (|3.15[l 
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analyticity to the domains below 



a{z)= < 



ap{z) 



ap{z) 



{z - ap)^> 

—2iiTUpS{3i{ap—z)) 



z e {z : Up <^{z) < Op+i} n U 
zen+n Dap,e 

zen-n Dap,e 



(3.17) 



In much the same way, we split the formula for a into a holomorphic and a singular 
part: 

«t(^) = ^pH^){^-apf"'^~^^''e'""''Kp{z) , ze'H+nDap,e; 
aliz) = ap{z){z-apy^''^x-^^-e''^''^Kp{z) , zeU^f^Da,,,. 



There, 



Kp{z) = 



{z - ak) 

^^.ap ii _ a, + ii _ ^^y, ^i^J. ^ ^ _ ^y, ^ ' J 

(3.18) 

The analyticity of F {z) on U guarantees that Kp {z) is holomorphic on the disk -Dap,e- 
Lastly, we define i , the regularized version of d^/i around a^, according to 



(0 = at (^) - «P)"'^ > (0 = {z) {z - a 



So that, 



dp{z) . 



(3.19) 



(3.20) 



3.3 The first step x ^- $ 

Let $ be related to x by 

$ (z) = X (^) [« {z)Y 



(73 = 



1 

-1 



(3.21) 



Then $ [z) satisfies the following RHP: 

• $ is analytic in C \ M ; 

• V/c € 1 1 ; n] , there exists G GL2 (C) such that: 

^ = Mkiyl2 + g{z)B{z) + \z-ak\{g{z) + \)ol^ \ \ )}[a(2)] 
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• (z) G$ (z) = (z) , z € M ; 

The function (7 and the rank one matrix B (z) are as given in (I3.3p and (I3.2p . In 
particular g has a singular behavior at z — >• given by (j3.4p . a, the solution of the 
scalar RHP (I3.14p . introduces an additional singular behavior at z — )■ 0^: 

O (\z -akf^'"'''^) for z^ak,zen+ 
O |z — Ofcl ^ for z ^ Ofc , z € 7t_ 



Finally, the jump matrix for <I> reads 



and 



P(z) = [1 - ^7-1 (z)] (z) , (3.24) 
Q(z) = [^-1 (z) - 1] (z) . (3.25) 

Clearly the solution of the RHP for $ exists as it can be built out of x- Its 
uniqueness can be seen along the same lines as in proposition 3.1. Note that, because of 
the different possible analytic continuations of a to the upper /lower half-planes (I3.17p . 
P and Q will have different analytic continuations to the complex plane depending on 
the value of 5? (z). In particular, 

P(z) = af^(z)-— ^-^-y-— — [x(z-ap)]2^'' z € D,,,, n , 

^ (3.26) 

^2iTTUpr.{di{ap~ z))+ixap 

Q{z) = Kp{z) — a2(z) zeDap,en'H^. 

[x (z - ap)\ P 

3.4 The second step $ T 

We now perform a transformation on The resulting matrix T will have its 
jump matrices exponentially close to the identity matrix, except in the vicinities of 
the singularities of a. The jump matrix G$ can be factorized into a product of an 
upper by a lower triangular matrix: 



G$ = M^M^ . (3.27) 



The matrices (resp. M^) 



Mtiz) = (J Pi^W'^y (3.28) 
M,iz) = (3.29) 
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admit analytic continuations from the intervals ] — oo ; ai [, ] ai ; 02 [, . . . , ] a„ ; +00 [ to 
some interval depending domains in U DT-L^ (resp.C/ n These analytic continua- 

tions are different if one starts from different intervals. In the following, so as to avoid 
any confusion, for > 0, (z) should be understood as the analytic continuation 
of Af| (JR. (z)) from the interval containing (z). A similar statement holds for (z). 

We draw a new contour F = U F_ IJfc=i '^+^ ^ '^-^ U. It allows to defines a 
piecewise holomorphic matrix function T (z) according to Figl2l 



T = , 


' + 


T = ^M-c 


-p(n) 
+ 


, T = 


ai 








an 






T = $M^-^ 




, T = $M^-^ 


< 




T = $ 


F_ 


*r 



Figure 2: Matrix T and its associated contour F = F^- U F_ Q F^^ U F^''. 

k=l 



As readily checked, T (z) is continuous across K\IJ^]^ and hence holomorphic 
in the interior of this new contour. By construction, T has cuts on the exterior contour 
r_|_ UF_. The additional cuts along IJ^^^ T^±^ are due to the different analytic contin- 
uations for P and Q (and hence and M^) to the strips {z : < (z) < Ofe+i}^Zj 
cf ()3.26p . The matrix T solves the following RHP: 

• T is analytic in C \ F ; 

• V/c G 1 1 ; n] , there exists G GL2 (C) such that: 



T{z) = Mk { h + g (z) B (z) + \z - ak\{g{z) + 1)0 



1 1 
1 1 



T 











T+{z)M^{z) = 


T_ 




z G F+ 


T+ {z)Ml\z) = 


= T_ 




z G F_ 


T+ {z) TVC) {z) = 


T_ 




ZGF? 


T+ {z)N^'^{z) = 




(^); 


.Gr« 



I e [1; n] 



M{z) 
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There, the rank one matrix B (z) ()3.2p and the function g (13. 3p are as in the RHP for 
X- We remind that g has a singular behavior at given by (|3.4p . The matrix M is 
expressed in terms of a and M^/^ according to: 



M{z) 



' a^^Mt (z) z^(n+\ u^=irf ) n v 



af (z) z(i\n-\ u^=irL^M n c/ 



(3.30) 



It is readily checked that the matrix M (z) has a singular behavior at z ^ given 

by 



M(z + afc) 



O 



O 



. |5RK) |^[min{-K{i.fc),R{P^)) 



for z ^ , z G 
for z ^ ^ , z (^U- 



Finally, the jump matrices (z) are defined by 



(3.31) 



ivW (^ 



iv(') (z) 



1 n/ (z) e*' 
1 

1 



ni (z) e " 
and their entries read 

ni {z) 



hm M7^ (z - e) M+ (z + e) 

e-^0 ^ 
5R(e)>0 

= lim M_ (z - e) Ml^ (z + e) 

e-i>0 
!R(e)>0 



2 G r 



(0 



z G r 



(0 . 



[a; (z - ai 
Ki (z) e*^*^' 
e^^'"'i^/(z) 

[x(z - a^)]^ ' 



-2i7r!/ 



1) • 



(3.32) 
(3.33) 



The solution of the RHP for T clearly exists and its uniqueness follows from a 
similar reasoning to proposition 3.1. The matrices T and x ^^re thus in a one-to-one 
correspondence. 

Note that, apart from vicinities of the points aj, z G |1 ; n], the jump matrices 
for T are exponentially close to the identity. We have almost been able to recast the 
original RHP into one suited for the Deift-Zhou steepest descent [T9| : it only remains 
to build the parametrices around the Oj's. 



4 Construction of the Parametrices, last transformation 

We first construct the parametrix for the model RHP on a small disc Z)o,e of 
radius e > and centered at 0. This model parametrix will be the key ingredient of 
the parametrices around the Oj's. 
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4.1 The model parametrix. 

The model parametrix P is a solution to the following RHP: 

• P is analytic Do,e \ {r_|_ U r_} ; 

• P = I2 + O ( —] , dDo e uniformly; 

\exj 

P+ (z) N{z) = P. (z) , zeT+nDo,e 

p+{z)N{z) = p.{z) , zer.nDo,. ' 

The jump matrices of this model RHP read 

l25 ixz 



N{z) = \ ^ K{z) 

) 1 



N{z) = K{z) 



\ l2<5 

\xz\ 



1 



1) 1 



and the function K [z) is assumed to be holomorphic and non- vanishing on Do,e- Note 
that the boundary 9-Do,e of the disk Z)o,e is canonically oriented just as depicted in 
Fig. O The RHP for P admits many solutions. For instance having one solution P, 




Figure 3: Set of contours in the RHP for the model parametrix. 



one can build another one by multiplying P on the left by a holomorphic matrix on 
Z)o,e that is equal to I2 up to corrections that are uniformly an O {1/ex) on (?Do,e a-^d 
hence on the whole disk L'o.e- 

A solution to the above RHP can be built thanks to the following procedure. We 
first assume that K \s a constant. Then the Riemann-Hilbert problem for P can be 
solved by the standard procedure. One performs the transformation 



P(z) = G(C) 



62 



-0-3 



with Q 



xz 



(4.1) 
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so that the jump matrices for @ are piecewise constant, and is a solution of a 
RHP on Do^xf This RHP is solved explicitly in the limit xe +00 by the standard 
differential equation method |28) . It is then enough to go back to the original matrix 
P. Eventually, we get that 



P(z) 



- 6,1 + 2j;-ixz) ibu"^ {1 + J + 6,1 + 2j;ixz)\ L 

-i62i^'(l + 7-5, l + 27;-ix2;) ^' (7 + 5, 1 + 27; ixz) / {xzf"'''^' ' 

(4.2) 

Here '^{a,c;z) is Tricomi's confluent hypergeometric function (CHF) given in (IA.3p . 
We remind that it has a cut on M_. The piecewise constant matrix L depends on 6 
and 7, 

e^^e"*^'^^ -7r/2 < arg (z) < 7t/2 



. TT(S-f) / 1 

2 I ^_2i^5„i^^ I tt/2 < arg (z) < vr ^ 



L 



e* 2 Q ^ j -vr < arg (z) < -vr/2 

whereas the coefficients 612 and 621 also exhibit an additional dependence on K 

Lastly, one has det [P] = 1. 

Using the asymptotic behavior of Tricomi's CHF (IA.6|1 . one readily checks that 
P has indeed the correct asymptotic behavior. The jump conditions can be checked 
thanks to the monodromy properties of Tricomi's CHF. These are (IA.4P for the jump 
condition on and (jA.SP in the case of the jump condition on r_. Moreover, P has 
no jump across Z^o.e \ ] ~f j [ • the discontinuity in L is there to compensate the one 
of {zy~^'^''^ . Hence P is holomorphic on [— e;0]. The fact that det [P] = 1 can be 
seen as follows. We first assume that |5R(7)| < 1/4. Then, using the local behavior 
at z = of Tricomi's CHF (jA.lOp . we get that det [P] = o {z~^). On the other hand, 
writing the jump conditions for det [P], one can easily convince oneself that the latter 
function is holomorphic on C \ {0}. Its singularity at z = is thus of a removable 
type. As det [P] 1 when z ^ 00, we get that necessarily det [P] = 1. To reach 
the case of generic parameters 7, we fix z 7^ and invoque the fact that ^ (a, c; z) 
is a holomorphic function of a and c. It follows that det [P] (z) is holomorphic in 7. 
As it is constant in the region |3f?(7)| < 1/4, we get that it is constant on C. Hence, 
det [P] (z) = 1 for all z 7^ 0. Now we get that det [P] is bounded in every punctured 
neighborhood of z = 0. It thus follows that it cannot have any power-law singularity 
at z = 0, and det [P] = 1. 

In this way we have built a solution of the RHP for P in the case of constant func- 
tions K. In order to extend this solution to functions K that are holomorphic and 
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non- vanishing in some open neighborhood of -Do.ej it is enough to notice that replac- 
ing the constant K appearing in the formulae above by a holomorphic non-vanishing 
function K (z) on Z)o,e doesn't change the analyticity of P on Do,e \ (r+ U r_), nor 
its asymptotic behavior on the boundary (9-Do,e • As the jump conditions hold point- 
wise, these are also satisfied. Thence we get a solution to the general RHP for the 
parametrix. 

4.2 The parametrix around 

Let Pap be defined as 



^ hp - Sp-, -iCp) ib^^ (z) ^ (1 + 7p + 6p; <p) 

-ib^f^ (z) ^{1 + jp- 6p; -iCp) ^ (tp + Sp-, iCp) 



Lr. 



{Cpf"'"'-'"^^ 



(4.5) 

There we have set C,p = x{z — ap) and the second argument of the CHF's is implicitly 
assumed to be 1 -|- 27p. The piecewise constant matrix Lp reads 



^ <sp-,p) / 1 

^ ^ ' Q Q-2i-K5p-iiT^p 



. n{Sp+jp) ( g-2j7r(5p+j7r7p q 

e* 2 

' 1 



-7r/2 < arg (z — ap) < 7r/2 
7r/2 < arg {z — Op) < vr 

-TT < arg {z — ap) < — 7r/2 



(4.6) 



and the coefficients (z) and b)^^' (z) are 



6g (z) 



-nrfp 



r (1 - 7p + 6p) 



Kp{z) r(-7p-(5p 



b^S{z) = -iKp{z)e'^^^^^'^ '^^^ 



r (5p - 7p) 



(4.7) 



Proposition 4.1 The matrix Pap {z) plays the role of a parametrix around ap in the 
sense that: 

• T Pap is holomorphic inside of Dap^e , 

• p-i =l2 + (l/2;i-2«('5p)) uniformly on dDap,e- 



Proof — It follows from det \Pap\ = 1 that Pap is invertible and that 
p-i =* Comat {Pap). The fact that P"^ =l2 + (l/x^-^^^^f)) uniformly on dDap^e 
follows from the asymptotic behavior on the boundary of ODq ,, of the matrix P defined 
in the previous section together with the fact that the function Kp (z) appearing in 
the definition of 6^2'' i^) a-^id 621'' i.^) depends on x'^^p, cf (|3.18p . 
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As Pap and T have the same jump matrices on It^"* U r^^| Pi Da^^e, ^Pa^ 
be analytically continued to the punctured disk -Dap.e \ {o-p}- It remains to see that 
the singularity at 2: = is removable. 

The local behavior of T at z — )• Op as well as the one of CHF at the origin (c/ 
appendix (lA.lOP ) imply that TP^^ has at most a power-law singularity at Up. In 
particular, it cannot have an essential singularity at ap. The singularity can only be a 
pole of some finite order. To set aside this possibility, it is thus enough to check that 
TP^^ is bounded in the quadrant < arg {z — Up) < 7r/2. 

The addition formulae for the CHF (lA.OP allow to express the product PapM^^a~'^^ 

as 

^r(l-27p)^>-5p V r(l-5p-7p)e2->ci>(5p_^p,l_27p;iCp) y ^'^ 

Here, once again, the second argument of Tricomi's CHF is assumed to be 1 + 27p. 
Recall that the $ functions are regular at = whereas the ^ functions have a 

power-law singularity of the type O ^("p "^^^^ . Since {z) /Qp' = O (1) for z ap, 
we get that 



P„^Mf 1 [a^]-^-' = {O (1) + 0[\z- a^r'"'^) + O (log \z - ap\)} O ( | | ) 



(4.9) 



(4.11) 



Moreover, using the representation (j4.8p . one obtains that 

Pap (z) Mf ^ (z) a-'^' {zY Comat {B (z)) = O ( } } ) . (4.10) 
Hence, using the local estimates for T around Op, we get that 

Pap (z) (z) = {O (1) + 0{\z- apj-^^) + O (log \z - ap|)} O | J 

Here we have used that det [T] = 1 what implies that, for 2; — >■ Op, 

T~i (z) = M-^a-'^' |/2 + g {z)' Comat {B {z)) + \z - ap\ {g {z) + I) O \ J ^ | M^^ , 

(412) 

for some Mp G GL(2,C). Therefore, as K (27p) < 1, we see that Pap{z)T ^ {z) 
cannot have a pole at z = ap. The singularity at 2; = Op is hence removable and 
Pap {z) {z) is analytic on Dap^e- As det [Pap (z) (z)] = 1, we have that 

T (z) P-/ {z) =* Comat (P,^ {z) T'^ {z)) (4.13) 
is also analytic on Dap,e- □ 
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4.3 The last transformation T — )• 17 

The matrix 



Q = 

satisfies the RHP 

• is analytic in C \ ; 



T 



z£D. 



Z G C \ U Da„e 



• n = i2 + o (i/z) 


, z - 


-> oo ; 






' n+ {z) (z) 


= 


(^) 


zer+ ; 




n+ (z) M-^ (z) 


= 0_ 


(^) 


z G r_ ; 


• < 


J7+(z)iV« (z) 


= 


(z) 






n+ (z)N^^'^ (z) 


= 


(^) 






. n+{z)Pa,{z) 


= 


(^) 


z G 



(4.14) 



The solution of this RHP for 0, exits and is unique as seen by already invoked argu- 
ments. The newly introduced contours are all depicted in FigHl 



Q = T 



Q = T 




n = T 




n = T 



Figure 4: Matrix and new contours Sq. 



The jump matrix vq for Q is uniformly I2+O (l/x^ in the {^n) and {^n) 
sense, ie there exists a constant c such that 



Here, we choose the following (R) and (M) matrix norms 



(4.15) 



\Ais) 



tr 



^t(s)A(s) |ds| , ||A (5)11^00(5.^) = max I (5)11^^(5.^) 



(4.16) 
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We also remind that p = 2max |3? < 1. By using the matrix integral equation 

k 

equivalent to the RHP for we see that — I2 uniformly. Moreover the first 

corrections to are uniformly O (l/x^^''). As a consequence, I2 is the unique solution 
of the RHP for il, up to the uniformly O (l/x^^'') corrections. 



5 Asymptotics of the resolvent 

Recall that the resolvent I — R oi I + V can be expressed in terms of /+ and /_, 
cf (I2.19p . Having solved the RHP for x perturbatively in x, we use this asymptotic 
solution in order to obtain the asymptotics to the leading order of the functions f±. 
The latter yield the leading asymptotic behavior of the resolvent. 

5.1 The zeroth order approximants to f±. 
Definition 5.1 Let o-nd x^'p'^ be the matrices 

x'Hz) = M-^{z)[a^{z)r^ , 
X'riz) = Pa^{z)M-'{z)[a^{z)]-^^ . 

Let be the solution of the RHP defined in subsection 14.31 and where all the 
circles in the contour So^, as depicted on figlH have a radius e. Then, if x is the 
solution of the RHP defined in subsection 13.11 we have 



^^Hz)x{z) = x^'^iz) , z eUn{z -.^(z) e]ap + e/2;ap+i-e/2[}{5.1) 

^2eiz)xiz) = Xp'iz) ,zeDa^,2e- (5.2) 

In a sense that will become clear in the following, x^^ is the leading solution of the 
RHP for X when z is uniformly away from the singularities at the Op's and Xp"^ is the 
leading solution of the RHP for x when z belongs to the disk -Dap,2e- 

The advantage of using the solution of Riemann-Hilbert problems with two sizes 
of disks around the a^'s , i}2e and $1^/2 is that the solution x+ (z) defined by (|5.ip 
on M \ U^^]^ [ afc — e ; Ofc + e ] + iO+ and by (15. 2p on U^^^ [uk - e ; at + e] + iO+ has a 
smooth correction matrix $7 around the gluing point z = Qpzbe. This will simplify our 
forthcoming analysis when integrating the solution x around the points = ±e. If 
we would have used a single solution Q^, then we should have had derived additional 
estimates for the behavior of this matrix around = =be, as a priori it could exhibit 
a non-smooth behavior there. Thence, we circumvent additional complications. 

Proposition 5.1 Let 

( If (2) ^ _ ,» ( e+ (^) \ ( /fi (^) \-.tcc ( e+ u 
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Then 



^Ifej^j ) =Ke-)""^( "^'/^^ ) for zG{ap<^{z)<ap+i}nU , (5.4) 
with Up having already been defined in (|3.15p . Also 



f+,P (-2) ^ _ e'^^^-'^^-^^^x-^^"^ I e+ (z) Sj'^ {z) 



-1 







f%{z) ) r(l-27p)x> 1^ S?)(z)e_(z) 

for Z G -Da„;2e 



r (1 + (5p - 7p) ^> (-7p -5p,l- 27p; -<p) 
r (1 - 5p - 7p) $ ((5p - 7p, 1 - 27p; zCp) 



The functions cx£^ were defined in (I3.19P whereas ^ is Humbert's CHF (IA.7P and 
Cp = X (z - Op) . 

Proof — 

The proof of first formula (j5.4p is straightforward. We explain how to derive an 
expression for in a vicinity of ap. In the intermediary calculations we suppose 
that z belongs to the quadrant {z : ap < 3fJ (z) < ap + e} n The final result is 

however valid on the whole disk as can be seen through a direct computation carried 
out on the other quadrants. In order to lighten the notations, the second argument of 
Tricomi's CHF is undercurrented to be 1 + 27p and we have set Q = x{z — ap). 

fr{z)\ _ ( ^,{^^-5p--iC) i6g)(z)^(l + 7p + 5p;iC) \ 



(^) J \ -ibf^ (z) ^ (1 + 7p - Sp-, -iO ^ (7p + Sp-, <) 

= CPe^f^-(e'i>C''''ai-(z)e_(z))"''(^p"|^^^) °) (5.5) 
^ ( ^i7p-5p;-iO <7S?^(l + 7p + -5p;<) W ^ 

(1 + 7p - -^p; -^C) * (7p + ^C) y V 1 y ■ 

(p)_ r(i + ^p-7p) (p)_ r(i-7p-jp) 

" r(-5p-7p) ' - r(5p-7p) ' ■ ^'-'^ 

And finally using the recombination formulas for CHF (lA.9p . the claim follows. □. 

Remark 5.1 The function f^, resp. f±^p, o,re good approximates to f± in their 
respective domains of validity. More precisely, for z € M \ U^^-^ [ap — e/2 ; Op + e/2] 

(t<t;i:<:j)^(x(^)-x"«)(-<:>)^(o»(.,-.)(;i^g 

(5.7) 



With 
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and for z G Da^^e, 



(5.8) 



5.1.1 Uniform estimates for the resolvent 



Definition 5.2 Let i?o {S,, rj) he called the zeroth order resolvent. We define it in terms 
offf andf^^ as : 

Ro it V) = (0 - 1) ^~%^[^-~r^l^^ ^"^^ ' r?,eGR\ [j[a^-e-a,+e] 

-Ro ??) = (o- (0 - 1) 2m - r?) ' ^ + • 

Using the exphcit expression for and we get the local expressions for the 
diagonal zeroth order resolvent 

2^^^ =i^-^i (log "+"-) (0 , (5.9) 

n 

where ^ gM \ IJ [ap — e;ap + e] and we remind that a+ and a_ are the boundary 
p=i 

values of a from the upper/lower half plane. Whereas for ^ G ] Op — e ; Op + e [ 

2me-^-ap (0 ""'j'^^i^'/^ = -^^^ (Tp, -^p; Cp) - (log <^ (7p, Sp- Q • 

(5.10) 

One should keep in mind that (p = x{^ — ap) and that ci.^\. (resp. Sj^l) is the 

boundary value of (resp. S|^^) from the upper (resp. lower) half-plane. In the 
above formulae we have introduced two functions 



r (7, S; t) = (-7 -5,1- 27; -it) $ ((5 - 7, 1 - 27; zt) 



1 - 27, 1 - 27 
1 + 5 — 7, 1 — (5 — 7 

+ (a,^>) (-7 - 5, 1 - 27; -it) $ (5 - 7, 1 - 27; it) 

+ $ (-7 - 5, 1 - 27; -it) (9,$) ((5 - 7, 1 - 27; ^^) , (5.11) 

and 

1 + 11^-7)^ = ^ (-'^ - 1 - -^0 $ (5 - 7, 1 - 27; ^t) (5.12) 

Also, we have used the standard notation 



nr(a,) 

oi,---,an \ fc=i 



61,..., 6„ 



m 

nr(6,) 

ifc=l 



(5.13) 
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We now focus on the relationship between the exact resolvent R and the zeroth 
order one Rq. Observe that one can write the exact resolvent R as 



R it V) = i-e- (0 (0) [x?^ (e)] " (,) ( I- ) , (5.14) 

Where x^^^ iO = ^^'^ (0 X (0 the leading order solution of the Riemann-Hilbert 
problem corresponding to the choice of the radius e' for the disks dDa^^e' centered at 
the singularities of the symbol ct(^). We remind that the subscript e in $7 indicates 
that the matrix $7 is the solution of the RHP corresponding to the contour of 



section 14^31 where all the circles centered at Op have a radius e. 

Note that the exact resolvent R rj) does not depend on the choice of e'. Hence, 
we can chose different values for e', depending on the point (^,7?) where we want to 
estimate the resolvent. According to the above remark, we can present the exact 
resolvent as 

R{i,v)=Ra{i.r^)+Rc{i,v) ■ (5.15) 

There we have introduced the correcting resolvent Re- This resolvent is defined in 
terms of matrices Vt^i where e' takes different values depending on the point where we 
are placed. More precisely. 



Rc {tv; 2) = {-f- (0 (o) ' — ( ;;;; i (5.16) 



for , ^ G R \ |J [ Op — e ; Op + e ] , and 
p=i 



RA(,v;2.) = (-/- K) /-(?)) "?.<0";-W-^^ ( /| W ) (5.17) 

for ?7 , ^ S ] Op — e ; Op + e [. The crucial point is that the correcting resolvent is indeed 
small in the sense of the integration: 

Proposition 5.2 Let fip G {^p-,lp}- Then 



Proof — It is a standard fact that a matrix 17^ (z) approaching the identity at oo and 
having a jump $7+:e {I2 + A) = on S^j, can be expressed in terms of solutions of 
a singular integral equation: 

iO = h + C^n. [n+-eA] (e) where n+., (C) = h + [57+;eA] (C) . (5.19) 
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It is readily seen from the asymptotics of CHF (lA.SP that the jump matrix I2 + A for 
J7e is such that ||A (-2)11^00(2^ ) = O (x''"-^), but also | |A| 1^^2(2^; ) = O (x''^^). In the 
above equation, we have introduced the Cauchy operator 

C^n. if] iO = I ^-fr^/ (-) > C^n. [/] (^) = C^n. [f] iO . (5.20) 

^ — )• means that ^ approaches € Sq^ non-tangentially from the left side of Sq^. 
It is a classical result that is a continuous operator on L2 (En J. The fact that 

II^IL2(Sq ) ~ 0[xf~^^, leads to | — -^2! |/^2(j]j^ ) = O [x^'^Y its turn, this 
means that the entries of C such that Qe{z) = I2 + C / z + O (-z"^) are a O (x^""*^) . 
We have that 



I i^.(C,C;./2)^ + ^^|^,(e,e^,,)^. (5.21) 

IR\ U [ap-e;ap+e] ""^ 
P=l 

There, we should substitute the values of the correcting resolvent Rc given in ()5.17p 
and (I5l6]l . 

The exact formula for the correcting resolvent is 

(e, e; e') = ^Y. (0 iO iO {n;,' iO d,n, (O),^ . (5.22) 

id 

There we agree upon fii = + and fj,2 = as well as on the fact that 

n 

f^HO eG^bk = K\ U [ak-e;ak + e] 



/f (0 = { 



k=l 



(5.23) 



fc=i 



Also e' = e/2 for ^ E Pbk and e' = 2e for ^ G 2^ioc- The choice of two possible values 
for e' depending whether ^ G T^hk or 2?ioc ensures that [Vt^i {^) d^Q^/ (C))jj is smooth 
on R. 

We first study the integral on M \ Up^i [up — e/2 ; Up + e/2]. As already pointed 
out, the jump matrix A is uniformly O (x''~^) . This ensures that Q {^) — 12 is also an 
O for ^ bounded. Moreover, by Lemma [6l 



1 1 / l + min(l,K) / I 1 



, for ^ ^ ±00 . (5.24) 
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As (0 ~ e^'^ for ^ ^ oo, and di3^a (^) = O (^"^) , we get that K22^ is absolutely 
integrable on M \ Up^^ [ Op — e/2 ; Op + e/2 ] and that the integral is a O . 

It remains to study the integral of (I5.22p on [ Op — e ; + e ] . ^l2e ) is smooth on 
this interval and equal to /2 + O Moreover, for ^ — >■ 



9/3,^^(6 = 0(le-apr'^Mog|e-ap 



(5.25) 



for /3p € {5p,7p}- The formula for Rc involves four terms. They can all be treated 
similarly so we only explain here how to estimate the contribution of the integral 
around Op involving the [/+;p] . We denote J+,+ this integral and have 



IJ4 



dp,^ (0 {^2' (0 d^n^. {0)21 \fZ] (0 de 



< ^p-2-2Ji{<5p)^/ 



-23?(7p) 



log(x/|z|)I«>(-5p-7p,l-27p;-iz)|2 dz . (5.26) 



We have used the bounds on all the smooth functions appearing in the first line and 
used the O estimates for dp^a (I5.25p . C is some constant that can depend on e and we 
have explicitly extracted the factor x^^^ from the bound on the product of Cl matrices. 

The integral in the last line of (|5.26p is divergent for x — >■ +00 so that its lead- 
ing X +00 asymptotics are obtained by substituting the z 00 behavior of the 
integrand. Since 

^ (Sp - 7p, 1 - 27p; -iz) = c±z^p+^p (1 + (1)) for z ^ ±00 , c± € C , (5.27) 
we have 



I |z|-2»(>hog(x/|z|)|cI>(-,5p-7p,l-27p;-^^)pd^ = y log (^) \t 



o 



[ex 



,2Ji(5p) + l 



loge 



1 



l + mm(l,^^) III 

1 1 



, ^ ±00 



(5.28) 



2K((5p) + l V '"''^ ' 2K((5p) + l 
Hence we have that |J++| < cx^"^ for some constant c, uniformly in 6p,^p. □ 
Lemma 6 Ve' > 0, 



(5.29) 



Here, we remind that k is such that F {^) = 1 + ( |^| 



when ^ — > ±00. 
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Proof — As already discussed, we have that Ve' > 0, (^) = I2 + o(l) in the 
L'^r\L°° sense. Hence, it remains to study the asymptotic behavior of d^^l^/ (^). 

We focus on the case ^ +00 as the ^ — )■ — 00 case can be treated similarly. We 
decompose the contour Eq^, = ^ U ^ , with 

S^^, =Sn^,n|z : K(z)<|| , Eg^ =Sf,. n|z : K(z)>|| . (5.30) 
Then, using the Cauchy integral representation for 0^/ (I5.19p . we get 

Where / + A stands for the jump matrix for Q^' and 

tR - 



ds 

^n+.^>A{s) . (5.32) 



2i7r {£, — s 

Reminding that ||^| 1^2(2 ,) ~ I \'^^\ i^) ^ ('^)] ' 



But taking the explicit formula for A (s) and using the estimates for the asymptotic 
behavior of 0^^,^^, (^), as well as the assumptions on the asymptotic behavior of the 

function F (^) = 1 + O ( |^[ for some k > 0, we get that 



tr{At(.)A(.)} = -^^— ^(l + o(l)) at .^00. (5.34) 

Where b = dist > uniformly in j^| large enough. Therefore, for such ^ 

+00 ^ ^ 
r \ds\ tr(At(s)A(g)| < 2 [ ds —-^-^(5.35) 



yR C 



e-2^^0 ( , ) . (5.36) 

1 ^l+mm(l,K) y ^ ' 



It follows that 



^l+mm(l,K) / ' 



for some constant c. The Lemma follows once upon replacing ||A||^2^2 , ^ ~ ^ ^'^'^ 
in dOT]) . □ 
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6 Asymptotics of the truncated Wiener-Hopf determinant 

In this Section we shall compute the leading asymptotics of determinants of trun- 
cated Wiener-Hopf operators that are Hilbert-Schmidt. 
Note that, for ^ G M, we can decompose (Tj^^^u^ — a^) into 

^v„F.(e-afc)=e^^(«)e^^(«). (6.1) 

Here, 

Qk (C) = 4 [log ii-ak-i)- log {i - Uk + + 2otS (a^ - 0] , (6-2) 

hk iO = Ik [log - flfc + + log - a/c - - 2 log - flfc |] , (6.3) 
and we remind that S is Heaviside's step function. It follows, 

dp^^a {z) = a {z) dp^ {gu {z) + hk (z)) for Pk G {Sk,lk} ■ (6.4) 

6.1 Integration 

Before carrying out the integrals appearing in the formulas for ds^ log det2 [/ + V2] 
and d^-^ logdet2 [/ + V2], we first establish a useful integration Lemma. 

Lemma 7 Let TZ be a Riemann integrable function on M, g( G 'la'^ {I,C), and I an 

o 

interval such that I D [ — ei ; €2 ], where ei > , £2 > 0. Then 

lim X€2 

j xdtg{t)n{xt) = g{0) j 7^(^) + o(l). (6.5) 



-ei — lim xe\ 



The 0(1) corresponds to terms vanishing in the ordered limit x — > +00 and then 

Proof — 
One has 

£2 £2 xe2 xy 

{g{t)-g{Q))n{xt)xdt = j dyg' {y) j dtKit) - J dyg' {y) j d^7^(^) (6.6) 

— ei xy —t\ —xei 

b 

Since TZ is Riemann integrable, we have that (a, b) ^ J IZ (t) dt is continuous and has 

a 

a finite limit at 00. Hence it is bounded, say by M > 0. It follows from (j6.6p that 

xdt{g{t) - g{0))n{t) < {ei + €2) M max l^'l □. (6.7) 

[-61 ;e2] 

We now establish the first separation identity. 
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Proposition 6.1 Let I + V be the generalized sine kernel defined by a symbol a — 1 
satisfying the Lp' (M) assumptions. Then the below identity holds 

(2...n) / , .X (2...n) , 

a| (ai + i) a)^ [ai — i) 
-"^ (ai + iO+) af -"^ (ai - iO+) 

V G(l -7i)G(l -7i) y 

There, G stands for the Barnes' function and the o (1) vanishes in the x +oo limit. 
The o (1) is uniform in 5p and 'jp. 

The functions cty^^^ refer to the values in the upper/lower half planes of the 
solution to the scalar RHP for a with vanishing exponents vi andVi, ie 



fc=2 ^ ft I / 



Finally, the notation \sj^=o means that we ought to set 5i = without altering all the 
other parameters. 

Proof — The asymptotic formulae for the resolvent (IS.lOp . (j5.9p combined with (j2.24p 
lead to 



ds, logdets [I + V] = I ||^^a5,a(0-95,F(e,0}d^ 



+ J 5,,a(0de. 

We have already shown in Proposition 5.2 that the second line is a O or more 

precisely, a term going to in the ordered limit x +oo and then e — > 0. We cannot 
send e to zero first, as the constant in the estimate for the O also depends on e, and 
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goes to infinity when e — t- 0. Writing the result of integration expHcitly yields 



ds^ logdet2 [I + V] = x J ^ {ds^ log a - ds^a} - J ^ {d^ log a+a_) (^) ds^gi (0 



ap+e 



^ I <j -xds^a (t) + ds^gi {t) 



2ix5r, 



p=l 



dt 



X {t — Up) + sgn {t — Up 
ixT (7p, 5p; xt) - f 9^ log aj^|a|^M {t + ap) (7^, 6p; xt) 



dt 
2^ 



+ / dt 



-IXT 



(7i,5i;xt) - (^(9^ log S|^|3|^l^ (t + ai)(^ (71,^1; + O [x^ ^) 

(6.10) 

There = M \ (Jp=i [op — e ; Op + e]. Also, we have introduced gi, the smooth part 
of 51, as well as r and the Riemann integrable regularizations of the functions r 
and ip introduced in (j5.12p . (jS.lip : 



51 (0 = f^i (log(C - ai - i) - log(C - ai 
r{^p,6p-xt) = e-*-'5-^s'^W|tr2>T(7p,5p;xt) + l 

^(7i,5i;xt) 



(6.11) 
(6.12) 



xt + sgn (t) 

^-i.Spssnit) \t\-^i.^(^^^^Sp;xt) - 1 (6.13) 



The Riemann integrability of r and if is part of the conclusions of Corollary IB. II given 
in Appendix ([B]) . 

The integrals in the third and fourth line of (I6.13P can be estimated using Lemma 
7. One gets that 



dt 



ixT {'yi,5i;xt) - f 9^ log a|^|s|"^M (t + oi) (71, 5i; xt) 



= -i J T{ji,6i;t)dt-x-^ (d^loga^^lafl^ (ai) j ^{ji,6i;t)dt + o{l) 
—00 —00 
= -26i + (71 + 61) V' (-71 - Si) + {Si - 71) ^ {5i - 71) + OT71 + o (1) (6.14) 

There we have used the value of the integrals of r (7, 6; t) and ip (7, 6; t) given in 
Corollary IB. II The o(l) stands for terms that vanish in the ordered limit x — > +00 
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and then e — )• 0. Very similarly 



^ / ^ [^51 51 (t + Op)] -ixT (7p, 5p; xt) - {d). log a|^|a|^l) (t + ap) (7^, 5^; xt) 



^ (flp) X -2vr7p = ^7^^5,5 (ap) . (6.15) 

p=i p=i 



The equality 



I id5,giit)dtlog{a^^la^fUt))^ = 0{e) (6.16) 

holds as we deal with a finite sum of integrals of piecewise smooth functions over 
intervals of length 2e. In order to estimate the integral appearing in the second line 
of (|6.13p . we need to change its expression a little 

ap+e 

\^ f f) 2ix(5p dt 

^ J '''^^^^^x(t-ap)+sgn(t-ap)27r 

^ n ~ 

r 25ix dt f sr- Q f 91 {t + flp) - ffi (flp - t) \ 2i5pXtdt_ 

J Xt-12'K J ^ t J Xt + 12TT 



2i7r 



P=i 

= -25i log (xe + 1) + O (e) = -25i log xe + o (1) (6.17) 

There we have used that (gi (t + Op) — gi {ap — t)) is smooth and xt/ (xt + 1) < 1. 
This ensures that the corresponding integrals are a O (e). The o (1) term stands, once 
again, for terms vanishing in the ordered limit x — ?• +00 and then e — )• 0. 

We now explain how to treat the integral containing dz log (a+a_) (z). There, the 
limit e — 7> cannot be taken directly as the integrand has non-integrable singularities 
at the afc. We thus start by integrating by parts: 

-J ^ (^C log a+a_) (C) ds^gi (0 = J ^ loga+a_ (0 d\^(.gi (0 

1 " ^ 

~ 2i7r X] ^-^i^-^i ~ ^) [a+aJ\ {ap - e) - ds^gi {ap + e) log [a+a_] {ap + e) 
p=i 

— log (a+a_) (ai — e) . (6.18) 

There, we have explicitly separated the regular part ds^'g of ds^g from the one contain- 
ing a jump. The sum appearing in the second line can be computed up to O (e log e) 
terms by using the local behavior of a+a_ around Op: 

a+Q_ {ap T e) = a^I'^ (ap T e + ^0+) Sj^^ (ap ^ e - iO+) e^^"^^'"^^^) \e\-'^^^ (6.19) 
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Hence, up to O (eloge) terms, only the discontinuous part of Q+a_ contributes to the 



sum: 



1 " ^ 

— X] {[dSiOSi loga+a_] (a^ - e) - [ds^giloga+a^] {ap + e)} 



2i7r 
p=i 

n 

= -Yl ^P^&i95, {ap) + O (e log e) . (6.20) 
p=i 

The singular behavior of a+a_ around oi leads to 

— log [a+a_] (oi — e) = 25i log e — m^i — log aP''"'"''' (ai + iO^) aj^''"'"'' (ai — iO"'') 



(6.21) 

^t/'i- have been defined in (16. 9p . 

Up to o (1) corrections, it is now possible to replace the integral over appearing 
in the rhs of (|6.18p by one over R as the integrand has integrable singularities at 
the points at- The resulting integral over R can then be evaluated by computing the 
residues at ^ = ai ± i thanks to the fact that a^/^ is analytic in the upper/lower 
half-plane and goes to 1 when z — )• oo in 'H±. One gets: 

j ^ log (a+a_) (0 dl ,^51 (0 = log («i + i) "-'■""^ («i " ^)) + 2^1 log 2 . 

(6.22) 



At the end of the day we get. 



/ ^ a+a-) (6 95,51 (0 = - log f 



(2...n) . X (2...n) , s 
(2...n) , (2...n) / ^ 

a+ («! + ?j a_ \a\ — i) 



+ 25i log 2e - i7r7i - ^ c/i (a^) 7p + o (1) . (6.23) 
Putting all the different results together we get 
ds, logdet2 [I + V2]=xl^ {ds, log a- ds,<j} + 26i log 0^ - 26i 

R 

/ (2...n) . X (2...n) , x \ 
- ^2:^^)7 , (2...n), T + + Si) ^ (-71 - 5l) 

+ (5i - 71)^^(71 -'5i) + 0(1) . (6.24) 

In particular, the two log e contributions from (|6.23p and (|6.17p cancel each other. We 
stress that o (1) stands for vanishing terms in the ordered limit x +00 and e ^ 0. 
It remains to integrate this result from up to 5i. This is licit as the remainders o (1) 
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are uniform in |K < 1/2. The integral of the ^ functions yields Barnes' functions 
due to the formula ()A.14p . Finally, we get that 



log 



det2 [I + V] 

= Si log 



x/§ log 



a 



a 



(2...n) . , (2...n) . .x 

a| (oi + I) [ai — i) 

-"^ (ai + iO+) (ai - iO+) 



+ %,.o}+^?log0) 



+ log 

This means that, for all e > 0, the limit 



^ G(l-7i + (^i)g(l-7i-^i) 1 , ... ..... 



det2 [/ + V] 



lim <log , , 
x^+oo I ^ \ det2 [/ + V] 



x/§ log 



0"! 



(6.26) 

exists. It is given by the rhs of (I6.25p . where o (1) are e dependent terms that go to 
when e — )■ 0. As the Ihs of ()6.25p is e-independent the x — )■ +oo limit cannot depend 
on e, therefore the value of the constant can be computed by sending e — )■ 0. The 
claim then follows. □ 

Proposition 6.2 Under the assumptions of the previous proposition, the following 
identity holds 



• detail +V\^^^^, ■ 



log 



+7i log [2)^^^ 



27r 



log 



(ai + z)af-"^(ai -iO+) \ ^ / G (1 - 71) G (1 - 71) \ 



(2...n.) 

a|^*""^ (ai + zO+) aj^*""'^ (ai — i) 



V G(l-27i 



+o(l) 



(6.27) 



Proof — Following very analogous steps to the (5i-derivative, one shows that 
d^, logdets [/ + " ^/ ^ {^^' log^Ui=o - ^71^^161=0} 

R 

/d^ " 
— [d^ log (a+a_ )] (0 /ii (0 + X] ^^-yi ^1 
„ — o 



p=2 



dt 



(1) 



+y^T^i97i/ii (z + oi) -i2;T(7i,0;3;t) - (^9^ log a| i j (t + oi) 99 (71, 0; xt) +o(l) . 



(6.28) 
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More precisely, we have replaced the integration of the exact resolvent by one involving 
the approximate one Rq for the price of O corrections. Then we have applied 

the integration Lemma 7 to estimate asymptotically the integrals around {ak}^^2 that 

n 

involve the CHF. These estimates produced the sum ^ Jpdyj^hi (op) and some o(l) 

p=2 

corrections. These, as before, are vanishing in the ordered limit x — )• +00 and then 
e ^ 0. 

The integral around oi should be considered separately as it is a little different in 
respect to the already studied integrals. We obtain 

/ 2i7r^'^'^^ + «i) 

—xe 

= 27i(?/'(l-7i)-2'0(l-27i) + l)-27ilogx + o(l) . (6.29) 

During the estimation of the above integral, one finds that the contributions stemming 
from the regular part of d^^hi only produce o (1) corrections as it vanishes at t = ai. 
Also one gets that the integral of the irregular part (equal to log |^|) versus ip produces 
at most 0(logx/x) corrections. The remaining integral in the second line can be 
estimated thanks to corollarv lB.il One should however use that fact that as t (71, 0; t) 
decreases at infinity as an oscillating factor dumped by 

ex 

f dt ^ , , , f dt ^ , , ^ / log x\ , ^ 

logx J — r (71,0; t) = log X J (^i, 0; + O f ^ j . (6.30) 

-ex R 

It now remains to study the limiting value of the integral in the second line of 
(I6.28p . In the case of the 71-derivative, this integral should be handled with greater 
care. Indeed, 

f d^ 

- / ^ (^5 log Q+a_) (0 9^1/11 (6 = 

- ^' / i (?TT - + Fw) ('"^ + ') - 

R\l-e;e] 

p=2 

+ /log _ - , } § . (0.31) 



-ixT (7i,0;xt) - f0^1ogS| j_S| i j (t + ai) (/J (71, 0; xi) 



(ttp - e) - log (a-^,+a4,_)^^-"^ 5^1/11 (op + e) 



40 



We have integrated by parts and decomposed the result into the integration of the 
singular part (line 2), and the regular part involving a|^'"''"''Qj^''"'"'' (lines 3 and 4). 
The latter functions have already been defined in (j6.9p . We have 

^ ^ 2ot(5(0 (6.32) 



in the sense of distributions ((5 (^) stands for the Dirac mass at zero). We can thus 
drop this functions from the first integral appearing on the rhs of the equation above. 
In particular, this integral is finite even if, a priori, it involves terms log \(,\ / (^ ± iO~^) 
that are integrated at distance e from zero. The sum appearing in the second line 
of (|6.3ip is handled similarly to the case of the 5i-derivative, ie by separating the 
smooth/singular parts of a||^'"''"^ around Up and then neglecting all the O(eloge) 
contributions. Finally, one can send e to zero in the integral appearing in the last line 
of (|6.31|) . This produces some corrections that go to zero with e due to the integrability 
of the integrand. The resulting integral over M can then be computed by the residues 
at ^ = ibiO^ and ^ = iti, exactly as it was done in the proof of proposition 16.11 At 
the end of the day, 

— {d^ log a+a-) (S,) d^.hi (^ = -27i log 2 - ^ 7p9^j /ii (op) 

Je 

/J2...n) ( (2...„) ( ^io+)\ 

+ log . (6.33) 

\^af-"^(ai+iO+)af-"^(ai-Oj 

Hence, 



— {9^, loga,,^^^ - 9^iCTU^^o} + 271 log 



«f..n) ^ (2...n) ^„^_io+)^ 

+log I ^ J 1 +271 (V' (1 - 71) - 2V (1 - 271) + l)+o (1) 

a| (oi + zO+) a| (ai - 



(6.34) 



We now integrate ()6.34p with respect to 71. The operation preserves the o (1) symbols 
as they are uniform in 71. The ip functions are integrated thanks to ()A.14p . Once 
upon integration, sending first x to infinity and then e to zero settles the value of the 
constant term. □ 



6.2 Asymptotics of the Predholm determinant 

Theorem 6.1 Let I+K he a truncated Wiener-Hopf operator acting on the segment 
[0;x] and generated by the symbol a — 1 with 

n 

(y{i) = F {i) Jl (C - Op) , fli e ai < • • • < On , (6.35) 
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where 

• F is holomorphic and non-vanishing in some open neighborhood of the real axis 
such that F -1 e L"^ (R), and even F (^) - 1 = O (^I^T^^, for some k > 0; 

• 5i(7fc) < 1/4 and \^{6k)\ < 1/2 , VA; E |1; nj. 
Then the leading asymptotics of det2 [I + K] read: 

n 

" (gfc) F'J (gfc) ^ / (g^ - + if \ , 

l}^ (afc + i) - 44 V ("'^ - + ("'^ - «f ) / 

(6.36) 

VFe /lat^e defined 

Gl \a\ = exp L I ^ [log (a) + 1 - a (0] i , (6.37) 

F [F] = exp I y" e [log F] (0 F-^ [log F] (-^ | • (6.38) 
We also remind that p = 2 max |K((^jfc)|. 

Proof — The result follows by a recursive applications of propositions 16.11 and 16.21 
At the end of the recursion, one also needs to invoke the Aheizer-Kac formula for the 
2-determinant of the truncated Wiener-Hopf operator generated by the (R) symbol 
F — 1 so as to fix the constant F [F] and the F dependent part of G| [a] . □ 

The leading asymptotics of Wiener-Hopf operators generated by general Fisher- 
Hartwig symbols (I6.36P reproduces all the previously know results: V A; i/^ = , or 
V /c z7fc = proven in [9j, V A; 7^ = proven in |12| and also the case of a pure 
Fisher-Hartwig singularity n = 1 and F = 1 |4]. We refer the reader to |11) for a 
restatement of all the know results in a language very close to the one used in this 
article. One should only pay attention to the different definition of Uuf^^v,, between this 
article and the book [llj. Indeed u^, resp. 17^, differ by an overall minus sign with 
respect to the conventions of the latter book. 

However, our result disproves the continuous analog of the Fisher-Hartwig conjec- 
ture [8] in its broad generality. Although most of the factors between the formula and 
the conjecture coincide, the latter predicts the presence of 

niT — — jV" • f^'^") 



k<p 



J«fc - apf + 4^ (gfc - apf ^ 
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whereas we find 



n( 

k<p \ 



{a-k - Op + 2z) (ofc - Op) 



(ofc - ap + i) 




(ofc - Up - 2i) (ofc - Op) 



(ofc - Op - i) 



) 



(6.40) 



The difference comes from the presence oiVpVk in the second exponent instead oiVkVp. 
Of course in all the cases previously investigated, the difference between the conjecture 
and the present result was not appearing as either the factor was not present or Up 
and Up were related by a sign. 

To end this Section we would like to stress that it is not a problem to obtain the 
sub-leading asymptotics of truncated Wiener-Hopf with Fisher-Hartwig symbols by 
the so-called x-derivative method: 



The above is a straightforward generalization of the identity for the pure sine kernel 
given in | l6j . The x +00 asymptotics for xi can be obtained by solving perturba- 
tively the singular integral equation satisfied by VL. We do not present the calculations 
here as we are going to derive the sub-leading asymptotics for the Toeplitz determinant 
case investigated in the next Section. 

7 Toeplitz matrices with Fisher-Hartwig type symbols 

In this Section, we adapt the previous analysis of the generalized sine kernel. In 
this way we obtain, in the framework of Riemann-Hilbert problems, the asymptotic 
behavior of Toeplitz determinants generated by symbols a having Fisher-Hartwig sin- 
gularities. This approach is based on an observation made by Deift, Its and Zhou 
in |16] concerning the relationship between the Fredholm determinant of a sine ker- 
nel on a circle and a Toeplitz determinant. Our results reproduce those obtained 
by T.Ehrhardt in his thesis [2'(J| . Moreover, the Riemann-Hilbert approach allows to 
compute sub- leading asymptotics to any order in a quite systematic, although quickly 
cumbersome way. At the end of this Section we shall establish the first sub-leading 
asymptotics of Toeplitz matrices with Fisher-Hartwig singularities. We observe that 
these sub-leading asymptotics of Toeplitz determinants partly restore the indepen- 
dence on the choice of a Fisher-Hartwig type representant for the symbol a. Indeed 
the jumps of a are characterized by parameters 5^- A shift of any 5^ by an integer 
describes the same jump. The freedom of choice of a Fisher-Hartwig representant 
for fj is broken if one considers the leading asymptotics only. These correspond to 
the choice < 1/2. However, part of the sub-leading asymptotics (the so-called 

oscillating ones) we obtain correspond to shifts (5j 1— )■ (Jj + 1 , 5^ 1— >• 5j — 1 in the pa- 
rameters appearing in the leading asymptotics. These sub-leading asymptotics shed a 
light on the mechanism appearing in the asymptotics for ambiguous case type symbols 
that has been conjectured by Basor and Tracy [3j and proven recently by Deift, Its 
and Krasovsky [14]. Indeed, terms that were subdominant in the asymptotic series 




for 2; ^ 00 . 
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for a generic set of parameters become of the same order of magnitude as the leading 
asymptotics when some of the parameters 5p and 7p are set to these specific ambiguous 
values. The global structure of the sub-leading asymptotics seems to follow the scheme 
already pointed out in |33| . We formulate a conjecture on this global structure at the 
end of this Section. Our conjecture can be seen as a generalization of the Basor- Tracy 
conjecture: we believe that the full asymptotic series for detm [T [a]] results of a 1 — 6 
periodization of only a small part of the asymptotic series. 



7.1 The Riemann-Hilbert Problem 

Let us consider an integral operator acting on the unit circle with the kernel 



V {z,z') = y/a{z)-Wa{z')-l 



m f—— —in t — 

Z-2 z 2 - Z 2 2;' 2 



2i7r (z — z' 



where 



a{z) = h{z)J[[l 



k=l 



1 



(7.1) 



(7.2) 



There we assume that b is holomorphic and non- vanishing in a vicinity of and has 
zero winding number. One can actually characterize the singular behavior of a on the 
contour of integration more explicitly. Namely, 



z 

1 

ak 



jSkie-dk-TTsgnie-ek)) 



(2 - 2 cos (6* - 9k)) 



Ik ' 



A = Qi{e-ek) 
ak 



We have set, just as in the preceding sections, 

25fc =yk-Vk 27fc = Uk + h'k- 



-6*^ G ] -TT ; vr [ . 
(7.3) 

(7.4) 



Here, we assume that |5R((5fc)| < 1/2 and K (7^) < 1/2. 

The authors of |16| observed that for m G N the pure sine kernel on the unit circle 
is of finite rank. This property persists in the case of the integral operator under 
investigation as: 



V (z, z') = ,/^{z)^l^a{z')-lY, 

p=i 



m ri— 1— — I — —V 



2z7r 



(7.5) 



Hence we have that 

det.^[/ + y] = detm-i 



1) e^^C^-^) 



det 



m—l 



2tt 



— afe^^U*'^'^-^) 
27r 



(7.6) 
(7.7) 
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We used the subscript ^ in order to insist that the Ihs is the Fredholm determinant 
of an integral operator acting on the unit circle ^ whereas the rhs is the determinant 
of an (m — 1) x (m — 1) matrix. 

Hence, the asymptotics of Toeplitz matrices with symbols having jump and power- 
law singularities will follow from those of the Fredholm determinant of the integral 
operator defined in (jT.ip . The only significant difference between the kernel (|7.ip and 
the one considered in the preceding Sections is the interval on which they act. Most 
of the steps in the derivation of the asymptotics are very similar. We only insist on 
the most striking differences. 

7.2 Asymptoic solution of the Riemann-Hilbert Problem 

We consider the RHP for a piecewise analytic matrix x having a jump on the unit 
circle 

• X is analytic on C \ ^ ; 

• V/c G 1 1 ; n I , there exists G GL2 (C) such that 

X = Mk S^h + g (z) B (z) + \z - ak\ {g {z) + 1) O (^^^ } ) | , ^ ^ Ofc ; 

• X h; 

2—)- 00 

. x+{z)G{z) = x-{z) ; ze^. 

There, just as for the Wiener-Hopf case, the rank one matrices Bi^ read 

Bk = ( ~L C ) ■ (7-8) 



1 

The function g is also defined similarly 



It has a singular behavior at z = Ofc of the type 

f 0(l)+0f(z-afc)"27A for 7fe/0 ^ 
g{z) = < ^' V ' ) whenz^ak- (7.10) 

( 0(log(z-afc)) for 7^ = 

We finally precise that the unit circle 'lo is oriented canonically [ie the + side of the 
contour corresponds to the interior of the circle) and that the jump matrix G reads 

We now define a new matrix T according to 



45 



• T = X a'^^ , for z being in the exterior of r_ and the interior of r+ ; 

• T = X , for z between r_ and ; 

• T = X al^^M^ , for z between r+ and 
Here a is the solution of the scalar RHP 

a analytic on C \ a- = aa+ , z G \ U^^;^ {ak} a ^ 1 when z ^ oo . 

(7.12) 

This scalar RHP can be solved explicitly in terms of the canonical Wiener-Hopf factors 
of 6: b = b+G [b] 6_. One has a = on Do,i and a = aj, on C \ -Do,i) with 

a^{z) = b-'{z)G[b]-'l[(l-^] (7.13) 

k=l ^ 

aUz) = b^iz)fl{l-^y'"' . (7.14) 



The matrices M-^^^ defining T read 
,^ , , /I (1 - a-^) a-'^ (z) z"'\ , , f 1 



1 '^'^^^-'-y {a-^-l)al{z)z-^ Ij ■ 

(7.15) 

We stress that, just as for the Wiener-Hopf case, the matrix M.^/^ (re*^) should be 
understood as the analytic continuation of M^-/^ (c*^) from a small neighborhood of 
e*^ to the ray [ e'^ ; re'^ ] . One readily sees that T satisfies the RHP 



• T is analytic in C \ F ; 

• Vfc G [1 ; n] , there exists G GL2 (C), such that: 



T{z) = Mk!^l2 + 9{z)B{z) + \z-ak\i9{z) + l)0(^'[ J ) | M (z) z ^ Uk 



T+ (z) Mt (z) = T_ (z) ; zG T+ 

r+iz)M-'iz) = T_{z); z€T_ ' 

rT,(.)^(0(.)^T_(.); .Gr|) ^ 

\ T+ {z) N^^' {z) = T_ (z) ; z G 

The function g and the rank one matrices B (z) are as defined above. Moreover, the 
matrix M reads 

a^'M^ z G |r»o,i \ ^k=i^+^} n U 

«rV ^e{c\(z)o,iU^=irf)}nc/ 
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The local behavior at z ^ ak of M can be inferred from the one of a^/^ and the 
explicit formulae for M-^/^. The jump matrices Ar(0 (2;), N^^^ {z) are defined as 



1 



and their entries read 



= _^(e2--*_i) . (7.18) 

[—imlog{z/ai)] ' 

We have defined, analogously to the Wiener-Hopf case, 

(7.19) 
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It is not a problem to see that the parametrix around Op can be chosen as 



* dp - ^p-, -Kp) 



* (tp + (^p; Kp) 



Ln 



-7p 



(7.20) 

Where we have set = —im\og{z/ap) and the second parameter of the CHF's is 
imphcitly assumed to be 1 + 27^. The piecewise constant matrix Lp reads 



e 2 e 



. 7r(^p-7p) 

e 2 



1 

Q g— 2i7ri5p— i7r7p 



-7r/2 < arg (Cp) < 7r/2 
7r/2 < arg (Cp) < tt 



(7.21) 







-TT < arg (Cp) < -7r/2 



and the coefficients {z) and 63^'' (z) are given by 



- Kp (.) r (-7p - 5p) 

The matrix 

n 

satisfies the RHP 



6g)(z) = -ii^p(z)e->^ 



r (1 - 7p - ^p) 

- 7p) 



p=l 



(7.22) 
(7.23) 



O is analytic mC\^n , Sq = (T \ 7) (J , with 7 = m U ^a^. 



,p=i 



• Q = 72 + (1/z) 


; z - 


-> 00 , 






' n+ (z) {z) 


= 




zeT+; 




n+ (z) M-^ (z) 


= n_ 


(z) , 


zeT_; 


• < 


n+ (z)iv(') (z) 


= n_ 


(z) , 






n+{z)N^^'^ (z) 


= 


(z) , 


.Gf« = 






= 


(z) , 


z G 



rj) \ <! r^^ n d, 



(0 



(0 



The solution of the RHP for 0,, clearly exits and is unique. Moreover it is uniformly 
l2 + 0{l/x^-p) with p = 2max|3i((5fe)|. 



7.3 The resolvent and asymptotics of logdet [/ + V] 



Just as in the case of kernels acting on M, one can reconstruct the approximate 
resolvent Rq in terms of Humbert's CHFs: 

1 — (z) ( m 
2m I z 



Ro {z, z) 



{-^ - dz^og{a+a-) (z)! 



(7.24) 
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forzG^\(U^n Da,,e 1 and 



\k=l 

Rq {z, z) = e-»^^p ^^'^ "+ ^^"^ i-—T (7p, 5p; -imlog z/up) 

-S^log (^Q^^S^^^ {z)ip{'yp,5p; -imlog z/ttp)^ +o(l) , 



n 

whenever 2; € |J ^ fl Daf.,e- The functions r (7, 5; z) and (7, 5; z) have been defined 
fc=i 

in (I51T]1 and (I5l2]l . 

One should observe that the proof of the differential identities (|2.23p is contour 
independent, so that we can also use them for a generalized sine kernel acting on 
Moreover one can prove that, just as for the Wiener-Hopf kernel, if R denotes the 
exact resolvent of the kernel (|7.ip . 

[ R{z,z) - Ro{z,z) 
dz ^— ^ dfs^a = o{l) . (7.25) 

The 0(1) terms vanish in the ordered limit x — t- +00 and then e — > 0^. Harping on 
the steps for the integration in the case of a kernel acting on M we get 

Proposition 7.1 Up to o(l) terms in the x — > +00 limit, one has 

^ f det[l + V] \ [ dz ( a \ 

, ^rr I T^i ] = m — — log - 5f log m 



«f '""^ (0) V . ^ G(l-7i + '5i)G(l-7i-^i. 1 , 
+ ^i^og , ^ ^ ,^ , ^ I +log — — — — +0(1) . 



4^;-") (ai) (ai) ) H G (1 - 71) G (1 - 71) 

(7.26) 



det[l + V\^^^^ \ f dz f CT\^^^, 



(2...n) 



/ 4,+ ("1) ^ 1 /G(l-7i)G(l-7i)\ 



(1) . (7,27) 



There, 



of-") (.) = 6- (.) G [6]- n (1 - f ) «f -"^ (-) = (-) n (1 - ^ 

k=2 ^ fc=2 



(7.28) 

^4/50, a^/| _|_/_ stand for their boundary values from the ± szrfes of^. 
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Proof — The proof goes exactly the same as in the case of Wiener-Hopf operators. 
The only notable difference is that one needs an additional version of the integration 
lemma. 

o 

Let e > 0, and / be a sub-interval of M such that V/ D [— e;e]. Then V/i G 
(/, C), and for any Riemann integrable on M function TZ (t) that it is decreasing as 
a power- law when \t\ oo, one has 

j mdtn (mt) h (t) log [2 (1 - cos t)] = h{0) j log [f] U (t) dt + o (1) . (7.29) 

-e R 

The o(l) term is vanishing in the ordered limit m +oo and then e — )• 0. The 
proof is straightforward by applying the original integration Lemma [7] to the function 
h (t) log [2 (1 - cos t)]-h (0) log [t^] that is (/, C). Moreover the power-law decrease 
of TZ (t) at infinity ensures that the contributions of the boundary which are of the 

+ CXD 

type logm J dtTZ{t) will indeed be subdominant with respect to o(l) . □ 



We point out that 



log = = / :^log ^ . (7.30) 



2i7rz \ cjk / J 2mz \ „ 

We chose to include these factors in the proposition above so as to make the parallel 
with the Wiener-Hopf case more obvious. By repeatedly applying Proposition 17. II and 
then invoking the strong Szego limit theorem we get 

Theorem 7.1 Let T[a\ he an m x m Toeplitz matrix generated by the symbol 

a(e) = b{z)ll(l--] (l-y) a,e'^,z = e'' (7.31) 

fc=l V ^ 

where b is analytic and non-vanishing in some open neighborhood of the unit circle ^ 
and 26k =^k ~ ^kj "^Ik = ^k ~\~ ^k ^'"e such that |5? (5^)1 < 1/2, (7^) < 1/2. 
Then the leading asymptotics of detm [T [a]] are given by 



detm [T [a]] = {G [b]) E [b] m^-^ \\ — — 

p=i ^ ~ 27p) 

n6r^(ap)fer'''(ap)n (i+^w) (^.32) 

r.-^ V "g/ 



where, 



2-K 

G[6] = e^^°s*l° E[b]=ek [log 6],,= / — e""' log 6 (61) . 

J 2tt 



(7.33) 

This reproduces the result of T.Ehrhardt |20| . One should pay attention that the 
exponents 5i and 7^ correspond to the exponents /3j and —Ui in Ehrliardt's notations. 
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7.4 The sub-leading asymptotics 

The authors of |16| found a way to express the discrete derivative of log [c] in 
terms of the RHP data. Their result reads 



Tm [c 



Xii (z = 0) 



(7.34) 



Where xii (0) stands for the upper diagonal entry of the solution to the RHP for x 
given at the beginning of this Section. Their proof also works, word for word, in the 
case of the more complicated kernel we consider, so we omit it here. It is now enough 
to determine the sub-leading asymptotics of the matrix Q defined by (|7.23p thanks to 
the singular integral equation: 



n (z) =h + 



ds 



2m {z — s) 



a+ (s) A(s) 



(7.35) 



with /2 + A being the jump matrix for il. The method for computing the corrections 
is standard. We send the reader to |.33j for more details. We stress that we did not 
chose this reference because of its originality with respect to the perturbation theory 
of such integral equations. We rather chose it as there, the perturbation theory is 
applied in notation quite close to the ones we use. 

It is easy to see that the jump matrix A is exponentially vanishing with respect to 
X away of the disks dDa^^e- However for s € dDap-e it admits the asymptotic expansion 



M 



1 



Af (s) + O {m-'^+P-^) 



^ ilm^ [log {s/ap)f 
We have set, using (a)„ = F (a + n) /T (a), 



(7.36) 



(s) 
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Up 



7| 



V 



§2 _ ^2 

ip 



I 



ilp - Sp)(, [--ip - 6p)^ 







-lY{-fp + 6p)^ (5p -7p)^ 



(7.37) 



and the functions {^) ^^'^ ^21 (^) defined in (|7.22p exhibit a slight dependence on 
m that is a O {mP). The standard considerations of a perturbative approach to ()7.35p 
lead to 



Vt (0) = h + 



Vli (0) , (0) , ^ 1 {a^/ajY'm^P 



+ 



m 



+ 



Where 



^^1 (z) 



^ A? jap) 
p=i 



z/ap ' 



(7.38) 
(7.39) 
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and the expression for 0,2 (0) is already more involved: 



d A 



(p) 



p=i 



ds 2s 



s — a-n 



log (s/ap) 



_^ ^1 (Qpj '^i (Q< 
^ l-a^/flp 

(p) ^ 

slog (s/ap) 



(7.40) 



We get 



Tm-l [a] 
This leads to 



/aj) m 



3p 



(7.41) 



logT^ [a] = mlogG H+logm V (7^ - (1 + o (1))+^ (1 + o (1)) . 

p=l 

(7.42) 

We recover the first two terms of the asymptotics of log [c] given in (|7.32p . The 
discrete m derivative method does not allow to determine the constant K, but this is 
irrelevant in what concerns the structure of corrections. Osc stands for what we call 
the oscillating corrections, whereas Nose for the non-oscillating ones. More explicitly, 



^ V -li + Si,-ip-5p J 



6+ (a^) b- (ap) {ap/ae - 1) {ae/ap - 1) 



n (ar/ae - l)-(>+'^'-) (a,/a, - l)^^-'" 



r=l 



n {ar/ap - l)-(>+'5'-) {ap/ar - ir^-'^ 

r=l 

T^P 



(7.43) 



whereas 



n 1 " 

Nose = ^ ap (7p2 - dl) [d, log ifp {z)l^^^ + 2 ^ ^'^^ " " + 

p=i p=i 

4t('^^^-7l)('^^7p^)27?^ (^-44) 

The oscillating corrections have a very nice relationship with the leading term 
T^^) {{ji} , {6i} ;m) [6] of the asymptotics given in (|7.32p . Indeed, one readily checks 
that, 



Ose = V t(0) (h} , {{6,}^^ ,6p + l,5,-l];m) 

p^g \ I p.? J / 



(7.45) 
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with b'pg (z) = a-pb (z) jaq. Such simultaneous changes (5p, Jg, 6) i-^ [bp + 1, (5^ — 1, 6pg) 
leave the value of the symbol a given in (j7.2p unchanged. Hence, this gives strong 
arguments supporting the Basor- Tracy conjecture as already a small part of the differ- 
ent Fisher-Hartwig representations for a is present in the sub-leading asymptotics for 
Tra \<y\ ■ However, our computations do not allow to give the proof of the Basor- Tracy 
conjecture in some particular cases where some terms in Osc become of the same order 
than T^^'^ ({7i} ; {^j} [^]- Although we formally reproduce some particular results 
of the conjecture (for instance = 1/2 = for some j and k) we cannot 

consider this limiting case as we do not have a control of the remainder. However, we 
raise the following generalization of the Basor- Tracy conjecture: 

Conjecture 7.1 The sub-leading asymptotics of a Toeplitz matrix with Fisher-Hartwig 
symbols slowly restore the broken by independence with respect to the choice of a 
Fisher-Hartwig representation for a . More precisely the asymptotics have the structure 

n 

^mM ~ J] ^({7^},R + ^i};"i)[fenJ , (z) = n M^) • (7-46) 

Sni=0 

With 

f ({7.} , {5.} ; m) [b] ~ ({7.} , {5.} ; m) [b] + g ^^^^^^^^^ (7.47) 

V k=l J 

having no oscillating terms with m. Note that the notation ~ stands for the equality 
in the sense of asymptotic series. 

This conjecture is a natural extension of the v periodicity conjecture raised in [33j, 
and, of course, of the Basor- Tracy conjecture. We stress that one could raise a similar 
type of conjecture in what concerns the Wiener-Hopf case. 



8 Conclusion 

In this article we have proven the formula for the leading asymptotic of Fred- 
holm determinants of truncated Wiener-Hopf operators generated by symbols hav- 
ing Fisher-Hartwig singularities. As a byproduct we reproduced, in the framework 
of Riemann-Hilbert problems, the leading asymptotics of Toeplitz matrices having 
Fisher-Hartwig singularities. We were also able to compute the first sub-leading 
asymptotics of Toeplitz matrices having Fisher-Hartwig singularities. These give sup- 
port to the Basor- Tracy conjecture. We proposed an extension of the latter conjecture. 
Our results were based on a connection between Toeplitz determinants and those of 
Wiener-Hopf operators: both are related to the so-called generalized sine kernel acting 
either on the unit circle ^ or the real axis M. In the case of Fisher-Hartwig singulari- 
ties, this generalized kernel has some jump discontinuities and power-law singularities 
on the contour. 
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An open question is the construction of an explicit asymptotic resolvent for trun- 
cated Wiener-Hopf operators I + K generated by symbols having Fisher-Hartwig type 
singularities. Indeed the resolvent is known in the Fourier space, so it would be enough 
to take the inverse Fourier transform so as to have the resolvent of / + i^T. Also, it 
would be interesting to find a proof for the conjecture we have raised above. 
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A Some properties of confluent hypergeometric function 

Tricomi's confluent hypergeometric function ^' (a, c; z) is one of the solutions to 
the differential equation [5]: 

zy" + (c-z)y'-ay = (A.l) 

For generic a and c, ^' (a, c; z) has a power-law singularity at the origin, and a cut on 
R~. It can be defined, for instance, by its Mellin-Barnes type integral representation 
in terms of Euler's F function 

7+ioo 

^'(a,c;z)= / F ' ' , ] , (A.2) 

^ ' ' ^ J V a,a-c+l ) 2m ' ^ ' 

'y—ioo 

that is valid for -3? (a) < 7 < min (0, 1 - 3? (c)) and -37r/2 < arg (z) < 37r/2. The 
latter integral representation is then supplemented by an analytic continuation. In 
the above formula we have used the standard hypergeometric type notation 

n 

. nr(afc) 

ai, . . . ,a„ \ _ fc=i ( \ n\ 

b^ b I ~ ■ ^ ' 

k=l 



Tricomi's CHF satisfies the monodromy properties 

line'" 
T {a,l + a — c) 



* (a, c; ze^-) = e'^-^M, (a, c; z) + — -— (c - a, c; -z) (A.4) 
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for < and 

^ (a, c; = e^-"^/ (a, c; z) - — -— (c - a, c; -z) (A.5) 

i (a, 1 + a — cj 

for Sz > . ^ (a, c; z) has an asymptotic expansion at z — ?> oo given by 

^ (a, C; = y (-1)" « («-C+l)n ^-a-n (^-M-aN ^ _ 3vr ^ ^^^^^^ ^ SvT ^ 

n! ^ ^ 2 2 

n=0 

(A.6) 

Humbert's CHF $ (a, 6; z) is another solution of (lA.ip . <I> (a, c; z) is an entire function 
that is defined in terms of its series expansion around z = 

+ 00 , . 

<P{a,c;z) = Y,^rrtz^, (A.7) 
It has the asymptotic expansion around oo 

a M 



r(c-a) V ^ / ^ nl{-z) V 



-a-M-l 



+ £(4e^z«-^ y "V^ + O (le^z'^-i-'^-^l) . (A.8) 

To ^ n!z" ^' M V ; 

^ ^ n=0 

There are many relations between these two different CHF. In particular 

^,(c-a+i)m^i-c^ (a - c + 1, 2 - c; z) = 

{r ( 1 r ( 1 !,) e-^*(c - ...c; -.)} (A^9) 

where e = sgn (Qz) and it is assumed that arg (z) G ] — 7r/2 ; 7r/2 [. One can also express 
^ (a, c; z) in terms of $ (a, c; z) thus allowing to access to the singularity structure of 
Humbert's CHF at the origin: 

^ {a,c; z) = T ^^'"^ ^ {a,c; z) + T (^"^ ~^ ^ z^-'<^ {a - c + 1,2 - c; z) . 

(A.IO) 

Actually a CHF is some limiting case of the Gauss hypergeometric function. This 
function is one of the solutions of the hypergeometric equation. We recall its series 
expansion around z = 0: 
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The above solution is regular around z = and can be continued to large value of z 
thanks to the identity 



a,b \ f c,b-a\, f a,l + a-c _i 



a,c — b J \ 1 + b — a 

One can also consider multi-variable generalizations of hypergeometric functions. We 
give here the definition of the Appell function of the second kind in terms of a double 
series that is convergent provided that \y\ + \z\ < 1: 

a- = > — — — y z . A.13 

We finally point out that the Barnes G-function admits an integral representation 
in terms of ip, the logarithmic derivative of Euler's Gamma function. 

^(z + l) = (27r)texp| - ^^^"^^ + y t^(t)dtl ^{z)>-l. (A.14) 



B Integrals of CHF 

B.l Series Expansion of Appell function of large arguments. 

Using the Mellin-Barnes type integral representation for CHFs [5j, Erdelyi was 
able to evaluate the Laplace transform of products of CHFs in terms of Lauricella's 
function [22]. The Lauricella function associated to an integral involving a product of 
two Humbert's CHF is better know as the Appell function of the second kind (jA.lSp . 
In terms of this function, Erdelyi's result reads: 

+ 00 

f dt e"*** 

f {p,ai,a2,€;s) = / ^27+i-e-p '^{S -lA - 27; -iait) ^ ((5 - 7, 1 - 27; ia2t) 


= s'^^-P-'T{e+p-2^) F2(p + e-2r, ~ 



-27 , 1-27 ' s' s ) ■ ^ ^ ^ 

Such integrals have been considered in |32) , in the case where e and p are integers and 
the integral is absolutely convergent. Here, we study the behavior of such integrals 
when e is close to zero. In that case, one cannot apply the integration procedure 
presented in [32] as it only applies to integer e. Moreover, we consider products of 
CHF that aren't decaying sufficiently fast at infinity. Hence one should regularize the 
integrals before taking the s — )• 0^ limit. Once a regularization is performed, this limit 
can be computed thanks to a series expansion of the second Appell function around 
00 
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We study / {p, ai,a2,e; s) as it is the generating function for all the integrals that 
appear in the evaluation of the trace of the resolvent. The precise procedure for 
computing these integrals will be explained in this Appendix. The idea is to derive a 
series expansion for F2 at 00 and then use it to compute, after a proper regularization, 
the s — )• 0"*" limit of / (p, qi, 02, e; s). 

Lemma 8 Let f (p, ai, Q2, e; s) be defined in terms of the Appell function of the second 
kind as in (jB.ip . Then f admits a series expansion around s = 0: 

— ^ — S 6 — ^ \ 

+ 52 (p, ai, 02, e; s) + S-^ {p, ai, 02, e; s) (B.2) 
Where the Si are series involving Gauss ' functions 

^, , e»f(p+e-27) fa^Y^ ( p + n-2j + e,p + n + e .s\ 

^n+p+e i-y z^^y^^j y + p + ^ + e 0,^7 

n + p + e-2^ 5 + j-p-e-n -5-j + n\ 
y 1 + n — 2j,l— p — n — €,l + n J 

s 



S2 {p, ai,a2, e; s) = . ^ 2-?^] 



n>l 



e''^ p / n + p + 5-7,7 + (5+p + n _._s_ 

gif 2^ 2i^M 1+p + n ' '02 



/ n+p + o — 7,n — e,e — — 7 — n 



(B.4) 



, s- \ 7+(5— p ■y—S 
7 + \ vra-i' 



S3(p,«.,a....)=e'-<'...e-i.r( ^7;_*^- ) ^^^^^^^^^^ 



n>0 



1 + 0— 7 — n,l + ny Vl + n- p — e 02 



n-7-(5 + e ^ $ f 5-7,7 + '^ ■ -i— 

l + (^-7-n-e,l + n + ey^"'^l 1 + n- p ' 02 



(B.5) 



Proof — We consider the Appell function of the second kind as in (IA.13|) and assume 
the following dependence between the parameters 

a = p + e + b + c and f = d=l + b + c = l + a — p — e (B.6) 

with e small and complex. The series expansion of the second Appell function can be 
re-summed into a Mellin-Barnes type integral representation [38j: 

a,6\„/ b c \ f ds f a + s,b + s,-s 

X 2Fi( ;y)(-xr . (B.7) 
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One possible choice of the contour ^ is depicted in fig. [6] and the integral is convergent 
provided |arg (— x)| + |$5 {y)\ < it. This can be seen using the asymptotics of F (z) in 
the region |arg(s)| < vr. One has 

r " + 5^6+^5 ' ) " ^'^^ ,a+f-/-ie-l^(«)l (1 + ,s^+ioo. 

(B.8) 

We remind the asymptotics behavior of a hypergeometric function of a large argument 



, a + s,c 
2Fi{ J. ■,y 



So that, all together 



^ { ^ ) {syy + 6^^+^^^ {syr-f^ [1 + {s-')] . 

(B.9) 



■ a + s,b + s,-s\ f a + s,c 



< c 



exp{-9(s)arg(-a;) + -7r|9(s)|} , (B.IO) 



where C is some computable constant. The integrand of (IB.7P is thus absolutely 
integrable provided |arg(— x)| + < tt. 




Figure 6: Contour of integration for the Appel function. 



We now split ()B.7p in two by using the analytic continuation of Gauss function for 
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for large y (fA?T2]l : 

' ds f a + s,b + s, —s,c — a — s \ f a + s,a + s + l — f l\ {—x) 

= I -r—l I , „ 2^1 



2i7r" V d + s, f — a — s J^'^K a + 1 + s — c ' y J 
. _ f <^s f a-c + s,b + s,-s,c\ ( c+l-/,c . -i\ i^-,.^ 

We were able to split the integral into two parts as each of them converges sepa- 
rately. The separate convergence of the integrals is readily seen from the asymptotic 
of the Gauss hypergeomertic function [5]: 



7 + s \ ^ + s 

-/3--(l + 0(s-i)) , \z\<l. (B.12) 



Hence putting 

, , ^ { a + s,b + s,c- a- s,-s \ ^ ( a + s,a- f + s + l 1 



d + sj-a-s J \ a-c + s + l ' y J {-yf ' 

(B.13) 

and using the asymptotic behavior of the Gamma function we get 

\g{x,y;s)\ < (7 e^'^l^('')l~^'^*K'^''s(-iO-arg(i-'(/-i)-arg(-|/)) ^ (B.14) 

Li is thus convergent in the region defined by the equation 

|arg (— x) — arg (l — — arg (— y)| < vr , and |?/| > 1 . (B.15) 

Similar calculation lead to the conclusion that L2 is convergent in the region 
|arg (— x)[ < vr. 

Li can be computed as a sum over the poles located at the right of ^. These are 
s = n , n G N and s = —b — e + n , n € N*. One eventually arrives to 

^ _ sin 7re )-'^ \^ ( Z^\"^ y ( a + n,b + n,p + e + n 

^ ~ ~sm7r(b + e)^~^' ^\y) \ d + n,l + n + e + p + b,l + n 

n>0 

( a + n,p + n + e .A sinvrft .-b-e , \-p-csr ( 

'""'[l + b + e + n+p-^y J- sin.(6 + e) ^""^ ^'^^ ^\~) 

\ / \ J n>l 



p + c + n,n - e,p + n - 6 \ ^ ( p + c + n,p + n - 6 ^ _i 



(B.16) 
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Similarly L2 can be computed as a sum over the poles located to the left of these 
are {—6 — n, —h — n — e} where n G N. This becomes apparent when one normalizes 
the Gauss function: 



«;^.)=^(c).$,(«;^.) , (b.it) 



so that 2^1 is an entire function of the parameters a, b and c. The result reads: 

{—iy{—y)~'^simT{d-b)-s-^{-xy f b + n,n — c\ ( —b,c 



\ ^ -\-n / ^ ^ 



(—2;)''+^ siuTre ^ \ 1 + n y yl + n — 

sin7r(d — 6 — e) / 6 + n + e, e + n — c\ / -5, c 



sin 77 {d — b) 



f^ + V''^""'! 2$if (B.I8) 



The joint condition for the convergence of L\ and L2 defines an open subset O of 
C^. Hence we can continue the series representation for the second Appell function 

to the largest open subset in containing O where the series is convergent. In 
particular, the series representation is well defined for the range of parameters that 
we use. Specifying the values of a, 6, c, d, f to the ones of the Lemma we obtain the 
claimed result. □ 



B.2 Useful integrals 

We will use the series expansion for F2 in order to compute some integrals of 
products of CHF. We remind the definitions of the functions r (7, 5; t) and if (7, 5; t) 

r ( 1 + } I f_ ^ ) ^ (7, t) = (-7 - 1 - 27; -it) $ (5 - 7, 1 - 27; it) 

+ {d,^) (-7 -5,1- 27; -it) $ ((5 - 7, 1 - 27; it) 

+ $ (-7 - (5, 1 - 27; -it) (5 - 7, 1 - 27; it) 

and 

r( ^l~y^'yj^_\^{i,5-t) = ^{-^-5,i-2r,-it)^{5-^,i-2r,it) . 



Corollary B.l Let \'R{5)\ < 1/2, Rij) < 1/2 and t {j,d;t), ip{-f,6;t) be as above, 
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then 




dt{T-2;z^{l.5;t)-l] =2i5 , dt { —-^—^ {^,6;t) - 1 ] = 



— oo 





( (7> + 1 - = -2i5-7r7+i (7 + <5) V (-7 " <5)+^ (5 - 7) (<5 - 7) 

r{-f,5;t) + l- J'^ I =-27r7 
i + sgn (t) y 

y" dtlog{\t\) (\t\-'^ ^{j,0,t) -1) = -27rj 
y Alog (itr^^ T (7, 0, t) + 1) = 27r7 (V (1 - 7) - (1 - 27) + 1) 

R 

The Riemann integrability of the integrands is part of the conclusion. 

Proof — Using the asymptotic behavior of Tricomi's CHF ()A.6P as weh as (IA.9P one 
can readily convince oneself that for < 1/2 

$ (-5 - 7, 1 - 27; -it) $ (5 - 7, 1 - 27;it) ~ 

t— >-+oo 



(5,$) i-5 - 7, 1 - 27; -it) $ (<5 - 7, 1 - 27; it) 

+ $(-5-7,1- 27; -it) {d,<^) (5-7,1- 27; it) 



t->+oo 



Where the terms that are sub-leading to e='=**t^2'^~^ are already absolutely integrable. 
For |3fi(5)| < 1/2, e^**t^^''~^ is only Riemann integrable, thus all integrals should be 
understood in this sense. Moreover the asymptotics in the region t — >■ —00 can be 
inferred from those at t — >■ +00 if one starts with the complex conjugated parameters 
6* and 7* and then takes the complex conjugate of the asymptotic series. This settles 
the question about the Riemann integrability of the different integrands. 

The proof of this Corollary is straightforward although quite long. The principle 
of the proof is to extract the divergent and constant terms in the s ^ 0^ limit from 
the expansion of / {p, ai,a2, e; s) around s ^ O"*". 

We shall explain in detail how to obtain the first four integrals. The remaining 
two are obtained in a similar fashion, although computations become more and more 
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involved. At the end of the proof we hst all the summation identities that are necessary 
to compute the s ^ 0^ limit in the other cases. Recall the notation introduced in 
(IB.ip . Then one has 

+00 +00 

y ^ (7, «, t) - t^V" } = lim^ J ^e- { V (7, 4, *) - t'-""''' } 

Since we compute the s — )• 0^ limit, it is enough to determine /(l,l,l,0;s) up to 
o (1) with respect to s — >■ 0"*'. One easily sees that Si (1, 1, 1, 0; s) = 0, c/ (IB.3p . Since 
we compute 52 (1, 1, 1, 0; s) (IB.4P up to O (s) terms, we can already replace Gauss' 
function by 1; the latter only contributes to higher orders terms in s. Then, for this 
particular choice of the constants p, ai,a2 and e, the second sum boils down to 

^2 (1, 1, 1,0; s) = -ie'^^ 5^rf^2)+0(s) = -^e*"^ + O (s) . (B.22) 

n>l ^ ^ 

Finally, we estimate 5*3 (1, 1, 1, 0; s) (jB.Sp . This term is the most complicated one. 
Indeed, due to the presence of the factor l/(ssin7re) in front of the sum, one must 
compute the linear in e terms of the sum (the zeroth order vanishes as it should be). 
Also, it is enough to expand the sum up to O {s log s) as such terms won't contribute 
to the result after the s ^ O"*" limit is performed. After some computations one gets 

53 (1, 1, 1, 0; = r ^ -V- 7 ) ~ + + ^ ■ 

(B.23) 

Adding up all the three contributions, we see that the part cancels with the one 
coming from the regularization term in (IB.2ip . The remaining terms combine to give 

+00 

j %[p (7, 5, t) - t27e-^} = -2z5e-^ . (B.24) 


Now, the first integral in the list of ()B.19P is obtained by considering (|B.24p in the 
case of parameters 6* and 7*, changing variables t — )> —t, and then taking the complex 
conjugate of the whole expression. 
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Wc now explain how to evaluate the integrals involving r {5, 7; t) . Since we have 
already established the value of integrals involving (6, 7; t), we only need to compute 



+00 



^ {(5,$) (-7 -6,1- 27; -it) $ (5 - 7, 1 - 27; it) 



+ $ (-7 - ,5, 1 - 27; -it) {d,<i>)i5 - 7, 1 - 27; it) 
2i5t^"<e^^^ ^( 1-27,1-27 



i + 1 Vl + '^-7,l-'^-7 
= -i lim |[5«J(0,l,l,0;s) -a„J(0,l,l,0;s)]e-^'^* 

In this case, only the first term appearing in Si contributes: 

- a«J .-Si] (0, 1, 1, 0; s) = ie--Tr ( ^ ~ J I ^ ^ ^ ) + O (s) . (B.26) 
Already for S2 one has to compute some less trivial sums 
[{da,-do,,).S2] (0,l,l,0;s) 

= ze-^rf^T^^'V'")lE ^ \ . \ J (i + Q(^)) 



—5 — 7, (5 — 7 y I n (n + (5 — 7) n (n + (5 + 7) 

.e-T ( 1 T V ) (^W-^d+^-l) ^ ^(l + ^ + 7)-V(l)| (1 ^ o 
\-5-7,()-7y|_ 7-(5 7 + (5 J' 

(B.27) 



And we have used >, — ^ = — — ^ -. Fmally, 

n>i n(n-a) a 



p„2-a„j.53] (o,i,i,o;s) = -e^-^r 



(5 — 7, —(5 — 7 
1 + 5 — 7, 1 — 5 — 7 

{2,5 [log {i/s) + 2i;{l)-i){-5--i)-il;{l + 5- 7)] + 1} (1 + O (s log s)) . (B.28) 
Adding together the three contributions we get 

e-^'^'r ( ^ +^-"27; 1 1 27" ^ ) [(^-^ - ^' ^' °' = 

- 2(5 [log (Vs) + V (1) - V' (-<5 - 7) - (1 + - 7)] - (7 + <5) V' ((^ - 7) 

2(5 TT ("(5 - e"*''(^+'') 

+ - + (7-(5)V'(l + .5 + 7)+ V ^ +0(glogg) 

7 — (5 sm7r(7 + ()) 

= 2(5 (log s - V' (1)) - i7r7 + (7 + (5) (-(5 - 7) + (5 - 7) - 7) + O (s log s) . 

(B.29) 
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We used the addition formulae for the ip function 

■0(1 + - V'(^) = ^ , 11^(1 + Z) -ll^i-z) = -TT cot TTZ , (B.30) 

in order to obtain the last line. The leading asymptotics of Tricomi's CHF around 
zero ^ (1, 1; s) = — log s + Tp (1) + {s log s) allows to take the s — )• 0^ limit in (IB.25p . 
We get 

+ 00 . 

« y ^|v?(7,'5;t) + T(7,<5;t)-^^ 

= — 7r7 — i (7 + (5) {—S — 7) — i (5 — 7) ((5 — 7) • (B.31) 







The other integrals involving r are then obtained from the latter results by the stan- 
dard manipulations that we have already described. 

The value of the last two integrals appearing in the Corollary is obtained by a 
similar procedure. Namely, 

^*^°^^{^(7,0;t)-t^n 



i27 







A long but straightforward computation yields 

r ( /J27! 1 I 27 ) (1, 1, 1, 0; = -7r7 + "^^^^"^"^^ + O log s) . (B.33) 

Similarly, 

+00 

dtlogt , , 

.5?+ -'^ ( 1 -"2^ L 27 ) [(^'^ - ^] °' ^) • (^-^^^ 







An even longer but as much straightforward computation leads to 

7 sin TT'j 



[{dl,-dlJS,]iO,l,l,0;s) = -^_^{27-V(l-27)-^(7) + i2} (B.35) 



[(Ce-c'aiO'^s] (0,l,l,0;s) = — (7rcot7r7-7-i) -^ + — -— cot2 7r7 

+ 7"^ [V- (1 - 7) + 27p (-7) - V (1)] - 27r7"^ cot 7r7 (1 + 27) 

H 7rcot7r7 — 7 (B.36) 

27 
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[{dl,e - da.e) (0, 1, 1, 0; s) = 7-' [ivr/2 + ^ (1) - ^ (1 - ^) - (-7)] • (B.37) 

Where we have dropped the O (s log s) symbol so as to lighten the formulae a little. 
In the intermediary computations of the contribution issued from S2 term we used the 
formulae below 



^ (7) - V' (-7) 



— 7^ 27 

n>l ' ' 



(B.38) 



4-^(n + 7)(n-7)' 27 

V^^^^ = -vr^{^(l + 27)+vrcotvr7} . (B.40) 
^ — ' — 7^ 7 

The last of these summation identities is maybe less standard. It follows from an e 
differentiation at e = of Dougall's formula for sums T functions [5]: 

f n + 7,n-7 + e \^ vr^ 

Y 1 + n — 7, 1 + n + 7 ) sin 7r7 sin vr (e - 7) r (1 - 27, 1 — 27 — e) ' 

(B.41) 

Finally, using standard properties of the -i/' function we get that 

- ( /_"27; L 27 ) - ^] 1' 0; ^) 

= 77r{V'(l-7)-2V'(l-27)} + 0(slogs) . (B.42) 
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